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ABSTRACT 

Using Gadget-2 cosmological hydrodynamic simulations including an 
observationally-constrained model for galactic outflows, we investigate how feedback 
from star formation distributes mass, metals, and energy on cosmic scales from 
z = 6 — > 0. We include instantaneous enrichment from Type II supernovae (SNe), as 
well as delayed enrichment from Type la SNe and stellar (AGB) mass loss, and we 
individually track carbon, oxygen, silicon, and iron using the latest yields. Following 
on the successes of the momentum-driven wind scalings (e.g. Oppenheimer & Dave 
2006), we improve our implementation by using an on-the-fly galaxy finder to derive 
wind properties based on host galaxy masses. By tracking wind particles in a suite of 
simulations, we find: (1) Wind material reaccretes onto a galaxy (usually the same one 
it left) on a recycling timescale that varies inversely with galaxy mass (e.g. < 1 Gyr 
for L* galaxies at z = 0). Hence metals driven into the IGM by galactic superwinds 
cannot be assumed to leave their galaxy forever. Wind material is typically recycled 
several times; the median number of ejections for a given wind particle is 3, so by 
z = the total mass ejected in winds exceeds 0.5f2{,. (2) The physical distance winds 
travel is fairly independent of redshift and galaxy mass (~ 60 — 100 physical kpc, 
with a mild increase to lower masses and redshifts). For sizable galaxies at later 
epochs, winds typically do not escape the galaxy halo, and rain back down in a halo 
fountain. High-z galaxies enrich a significantly larger comoving volume of the IGM, 
with metals migrating back into galaxies to lower z. (3) The stellar mass of the 
typical galaxy responsible for every form of feedback (mass, metal, & energy) grows 
by ~ 30 x between z = 6 — > 2, but only ~ 2 — 3x between z = 2 — > 0, and is around or 

1 /3 

below L* at all epochs. (4) The energy imparted into winds scales with M ^ , and is 
roughly near the supernova energy. Given radiative losses, energy from another source 
(such as photons from young stars) may be required to distribute cosmic metals as 
observed. (5) The production of all four metals tracked is globally dominated by 
Type II SNe at all epochs. However, intracluster gas iron content triples as a result 
of non-Type II sources, and the low-z IGM carbon content is boosted significantly by 
AGB feedback. This is mostly because gas is returned into the ISM to form one-third 
more stars by z = 0, appreciably enhancing cosmic star formation at z < 1. 

Key words: intergalactic medium, galaxies: abundances, galaxies: evolution, galax- 
ies: high- redshift, cosmology: theory, methods: numerical 



1 INTRODUCTION 

Galactic-scale feedback appears to play a central role in 
the evolution of galaxies and the intergalactic medium 
(IGM) over the history of the Universe. Mass feedback in 
the form of galactic outflow s curtails star formation (e.g. 
ISpringel fc Hernauisdl2003bl hereafter SH03b) by removing 
baryons from sites of star formation, thereby solving the 
overcooling problem where too many baryons condense into 



stars (e.g. Bal ogh et aU feoOll. The energy in these winds 
carry metal-enriched galactic interstellar medium (ISM) 
gas out to large distances, where the metals are observed 
in quasar-absorption l ine spectra tracing the IGM (e.g. 
ICowie fc Songailalll998l ). Galactic outflows appear to be the 
only viable method to enrich the IGM to the observed levels 
as simulations show t hat tidal stripping only is no t sufficient 
l|Aguirre et al.l l200ll ; lOppenheimer fc Pavel I200R hereafter 
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OD06). Hence understanding galactic outflows is a key re- 
quirement for developing a complete picture of how baryons 
in all cosmic phases evolve over time. 

Modeling galactic outflows in a cosmological context has 
now become possible thanks to increasingly sophisticated al- 
gorithms and improving computational power. The detailed 
physics in distributing the feedback energy from supernovae 
and massive stars to surrounding gas still remains far below 
the resolution limit in such simulations, so must be incor- 
porated heuristically. There are t wo varieties of ap proaches 
of feedback: thermal and kinetic. [Kobayasbj| (|2004h injects 
energy from galactic superwinds and supernovae thermally 
into a number of surrounding gas particles in a Smoothed 
Particle Hydrodynamic (SPH) simulation, and find that hy- 
pernovae with 10 x the typical supernova energy are needed 
to enrich the IGM to observed levels while matching the stel- 
lar b ar yonic content of the local Un iverse l|Kobavashi et all 
l200Tt ). ISpringel fc Hernquistl i|2003ah (hereafter SH03a) in- 
troduced kinetic feedback in SPH cosmological simulations 
where individual gas particles are given a velocity kick and 
their hydrodynamic forces are shut off for a period of 30 Myr 
or until they reach 1/10 the star formation density threshold. 
By converting all the energy from supernovae into kinetic 
outflows with constant velocity, SH03b are able to match 
the star f ormation history of th e Universe while enriching 
the IGM. ICen fc Ostrikerl | 20061 ) introduce a kinetic wind 
model in grid-based hydrodynamic simulations, and are able 
to match the obser ved IGM O VI lines in the local Universe 
jCen fc Fandl2006h . 

In OD06 we took the approach of scaling outflow prop- 
erties with galaxy properties, and explored a variety of 
wind models winds in Gadget-2 simulations. We found 
that the scaling s predicted by momentum-driven ga lactic 
superwinds fe.g. lMurrav. Quatert. fc Thomp son 2005;, here- 
after MQT05) provide the best fit to a variety of quasar 
absorption line observations in the IGM, while also repro- 
ducing the observed cosmic star formation history between 
z = 6 — ► 1.5. In the momentum-driven wind scenario, ra- 
diation pressure from UV photons generated by massive 
stars accelerates dust, which is collisionally coupled to the 
gas, thereby driving galactic-scale winds. MQT05 formu- 
lated the analytical dependence of momentum-driven winds 
on the velocity dispersion of a galaxy, er, deriving the re- 
lations for wind velocity, i> w ind oc a, and the mass loading 
factor (i.e. the mass loss rate in winds rel ative to t he star 
for mation rate), r\ oc a -1 . Observatio ns by iMartinl (|2005al ) 
and lRupke. Veilleux fc Sanders! |2005h indicate i; W md is pro- 
portional to circular velocity (where v c irc ~ a) over a wide 
range of galaxies ranging from dwarf starbursts to ULIRG's. 
Mass outflow rates are difficult to me asure owing to the mul- 
tiphase natur e of galactic outflows (|Strickland et all |2002| ; 
lMartin|[2005bl ). but at least at high-z there are suggestions 
that the mass outflow rate is of the order of the sta r forma- 
tion rate in Lyman break galaxies l|Erb et alj|2006l ). A the- 
oretical advantage of momentum-driven winds is that they 
do not have the same energy budget limitations as do su- 
pernova (SN) energy-driven winds, where the maximum is 
~ 10 51 ergs per SN, because the UV photon energy gen- 
erated over the main sequence lifetime of massive stars is 
~ 100 x greater l|Schaerer|[2003l ). OD06 found that trans- 
forming all SN energy into kinetic wind energy often is not 
enough to drive the required winds, particularly at lower 



redshifts. Moreover, galactic-scale simulations find that in 
practice only a small fraction of SN energy is transferred to 
galac t ic-sca l e winds |Mac Low fc Ferrar a 1999; Fuiit a et al.l 
12004 120081 ; ISpitoni et alj|2008l ). In short, the momentum- 
driven wind scenario seems to match observations of large- 
scale enrichment, is broadly consistent with available direct 
observations of outflows, and relieves some tension regarding 
wind energetics. 

Still, for the purposes of studying the cosmic metal 
distribution, the exact nature of the wind-driving mech- 
anism is not relevant; in our models, what is relevant is 
how the wind properties scale with properties of the host 
galaxy. The inverse a dependence of the mass loading fac- 
tor appears to be necessar y to sufficiently curtail star for- 
mation in high-z galaxies l|Dave. Finlator fc Oppenheimerl 
120061 : iFinlator et all 120071 ). At the same time, they enrich 
the IGM to the observed levels through moderate wind 
velocities that do not overheat the IGM (OD06). Contin- 
ual enrichment via momentum-driven wind scalings repro- 
duces the relative constancy of f2(Crv) from z ~ 6 — > 1.5 
(OD06) and the approximate amoun t of metals in the var- 
ious baryonic phases at all redshifts (iDave fc Oppenheimej 
120071 . hereafter DO07). The observed slope, amplitude, and 
scatter of the g alaxy mass-metallicity relation at z — 2 
llErb" et al.ll2006|) is reprod uced by momentum-driven wind 
scalings (jFinlator fc Davil200Sl ). While only a modest range 
of outflow models were explored in OD06, the success of a 
single set of outflow scalings for matching a broad range of 
observations is compelling. This suggest that simulations im- 
plementing these scalings approximately capture the correct 
cosmic distribution of metals. Hence such simulations can 
be employed to study an important question that has not 
previously been explored in cosmological simulations: How 
do outflows distribute mass, metals, and energy on cosmic 
scales? 

In this paper, we explore mass, metallicity, and energy 
feedback from star formation-driven galactic outflows over 
cosmic time. We use an improved version of the c osmological 
hydrodynamic code Gadget-2 (|Springell 120051 ) employing 
momentum-driven wind scaling relations, with two major 
improvements over what was used in OD06: (1) A more so- 
phisticated metallicity yield model tracking individual metal 
species from Type II SNe, Type la SNe, and AGB stars; and 
(2) An on-the-fly galaxy finder to derive momentum-driven 
wind parameters based directly on a galaxy properties. The 
OD06 simulations only tracked one metallicity variable from 
one source, Type II SNe, and used the local gravitational 
potential as a proxy for a in order to determine outflow 
parameters. These approximations turn out to be reason- 
able down to z ~ 2, but at lower redshifts they become in- 
creasingly inaccurate; this was the primary reason why most 
of our previous work focused on z > 1.5 IGM and galaxy 
properties. By low-z, Type la SNe and AGB stars con- 
tribu t e significantly to cosmic en richment (|Mannucci et al.l 
120051 : IWallerstein fc Knapp|[l998h . and these sources have 
yields that depend on m etallicity (jWooslev fc Weaverlll995l ; 
iLimongi fc Chi cffi 2005). Our new simulations account for 
these contributions. Next, using the gravitational potential 
wrongly estimates a especially at low-z, when galaxies more 
often live in groups and clusters and the locally computed 
potential does not reflect the galaxy properties alone (as 
assumed in MQT05). This tends to overestimate u w ind and 
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underestimate 77, resulting in unphysically large wind speeds 
and insufficient suppression of star formation at low red- 
shifts. Our new simulations identify individual galaxies dur- 
ing the simulation run, hence allowing wind properties to be 
derived in a manner more closely following MQT05. 

The paper progresses as follows. In §2 we describe in 
detail our modifications to Gadget-2, emphasizing the use 
of observables in determining our outflow prescription and 
metallicity modifications. §3 examines the energy balance 
from momentum-driven feedback between galaxies and the 
IGM using the new group finder-derived winds. We follow 
the metallicity budget over the history of the Universe in 
§4 first by source (§4.1), and then by location (§4.2), briefly 
comparing our simulations to observables including C IV in 
the IGM and the iron content of the intracluster medium 
(ICM). §5.1 examines the three forms of feedback (mass, 
metallicity, and energy) as a function of galaxy baryonic 
mass. We determine the typical galaxy mass dominating 
each type of feedback (§5.2). We then consider the cycle 
of material between galaxies and the IGM, introducing the 
key concept of wind recycling (§5.3) to differentiate between 
outflows that leave a galaxy reaching the IGM and halo foun- 
tains - winds that never leave a galactic halo. We examine 
wind recycling as a function of galaxy mass in §5 . 4. §6 s um- 
marizes our results. We use I Anders fc Grevessd i| 19891 ) for 
solar abundances throughout; although newer references ex- 
ist, these abundances are more easily comparable to previous 
works in the literature, and we leave the reader to scale the 
abundances to their favored values. 



2 SIMULATIONS 

We employ a modified version of the N-body+hydrodynamic 
code Gadget-2, which uses a tree-particle-mesh algo- 
rithm to compute gravitational forces on a set of par- 
ticles, and an entropy-co nserving formulation of SPH 
i|Springel fc Hernquisdl2002r ) to simulate pressure forces and 
shocks in the baryonic gaseous particles. This Lagrangian 
code is fully adaptive in space and time, allowing simula- 
tions with a large dynamic range necessary to study both 
high-density regions harboring galaxies and the low-density 
IGM. 

Gadget-2 also includes physical processes involved in 
the formation and evolution of galaxies. Star-forming gas 
particles have a subgrid recipe containing cold clouds em- 
bedded in a warm ionized medium to simulate the processes 
of evaporation and cond ensation seen in our own galaxy 
l|McKee fc Ostrikerl Il977l ). Feedback of mass, energy, and 
metals from Type II SNe are returned to a gas particle's 
warm ISM every timestep it satisfies the star formation den- 
sity threshold. In other words, gas particles that are eligi- 
ble for star formation undergo instantaneous self-enrichment 
from Type II SNe. The instantaneous recycling approxima- 
tion of Type II SNe energy to the warm ISM phase self- 
regulates star formation resulting in convergence in star for- 
mation rates when looking at higher resolutions (SH03a). 

Star formation below 10 Mq is decoupled from their 
high mass counterparts using a Monte Carlo algorithm that 
spawns star particles. In Gadget-2 a star particle is an 
adjustable fraction of the mass of a gas particle; we set this 
fraction to 1/2 meaning that each gas particle can spawn 



two star particles. The metallicity of a star particle remains 
fixed once formed; however, since Type II SNe enrichment is 
continuous while stars are formed stochastically, every star 
particle invariably has a non-zero metallicity. The total star 
formation rate is scaled to fit the disk-surf ace density-star 
formation rate observed by |Kennicutt] (fl998), where a single 
free parameter, the s tar formation timescale, is set to 2 Gyr 
for a ISalpeterl l| 19551 ) initial mass function (SH03a). 

Even with self-regulation via the subgrid 2-phase ISM, 
global star formation rates were found to be too high, mean- 
ing another form of star formation regulation is required. 
SH03b added galactic-scale feedback in the form of kinetic 
energy added to gas particles at a proportion relative to their 
star formation rates. They set the wind energy equal to the 
Type II SNe energy, thereby curtailing the star formation 
in order to broadly match the observed global cosmic star 
formation history. SH03b assumed a constant mass loading 
factor for the winds, which resulted in a constant wind ve- 
locity of 484 km/s emanating from all galaxies. OD06 found 
that scaling the velocities and mass-loading factors as pre- 
scribed by the momentum-driven wind model did a better 
job of enriching the high-z IGM as observed, while better 
matching the cosmic star formation history. 

We have performed a number of modifications to 
Gadget-2 since OD06. These include (1) the tracking of 
individual metal species, (2) metallicity-dependent super- 
nova yields, (3) energy and metallicity feedback from Type 
la SNe, (4) metallicity and mass feedback from AGB stars at 
delayed times, (5) a particle group finder to identify galaxies 
in situ with Gadget-2 runs so that wind properties can de- 
pend on their parent galaxies, and (6) a slightly modified im- 
plementation of momentum-driven winds. We describe each 
in turn in the upcoming subsections. 

All simulations used here are run with cosmolog- 
ical parameters consiste nt with the 3-year WMAP re- 
sults (Sp ergel et al.l 120071 ). The parameters are f2o = 0.30, 
n A = 0.70, D. b = 0.048, H = 69 km s _1 Mpc -1 , 
as = 0.83, and n — 0.95; we refer to this as the l- 
series. Note that og is somewhat higher than the WMAP3- 
favored value of 0.75, owing to observations suggesting that 
it ma y be as high as 0.9 fe.g. lRozo et alj|2007l ; lEvrard et alj 
2008). Our general naming convention, similar to OD06, 
is \(boxsize)n(particles/side)vzw- (suffix) where boxsize is in 
/i -1 Mpc and the suffix specifies how the winds are derived 
( "a" from the on-the-fly group finder, or from the local 
gravitational potential) and whether AGB feedback was not 
included ("nagb"). 

Table [T] lists parameters for our runs presented in this 
paper. We ran a series of test simulations with 2 x 128 3 
particles in 8 and 32 /i _1 Mpc boxes to explore the effect of 
turning off the AGB feedback and using the old prescrip- 
tion of using potential-derived i> W md- The 2 x 256 3 simula- 
tions are our high-resolution simulations and range in gas 
particle mass from 0.59 - 302 x 10 6 Af Q . The 132n256vzw-(T 
simulation was by far the most computationally expensive 
simulation taking in excess of 50,000 CPU hours on an SGI 
Altix machine. The 116n256vzw-cr simulation contains the 
minimum resolution needed to resolve C IV IGM absorbers 
(OD06), however it is prohibitively expensive to run this to 
z — 0. We will use the 18nl28vzw simulations at the same 
resolution but a smaller box to explore these absorbers; this 
box appears to converge with the U6n256vzw simulation at 
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Table 1. Simulation parameters 



Name" 






m SPH 


mi , 

dark 


M d 'f . 

gal, mm 




AGB Feedback? 


Wind Derivation 


Test Simulations 


18nl28vzw-<I> 


8 


1.25 


4.72 


29.5 


151 


0.0 


Y 


$ 


18nl28vzw-er 


8 


1.25 


4.72 


29.5 


151 


0.0 


Y 


cr 


18nl28vzw-cr-nagb 


8 


1.25 


4.72 


29.5 


151 


0.0 


N 


a 


132nl28vzw-$ 


32 


5.0 


302 


1890 


9660 


0.0 


Y 


* 


132nl28vzw-o- 


32 


5.0 


302 


1890 


9660 


0.0 


Y 


a 


132nl28vzw-er-nagb 


32 


5.0 


302 


1890 


9660 


0.0 


N 


a 








High-Rcsolut: 


ion Simulations 






18n256vzw-er 


8 


0.625 


0.590 


3.69 


18.8 


3.0 


Y 


IT 


116n256vzw-CT 


16 


1.25 


4.72 


29.5 


151 


1.5 


Y 


a 


132n256vzw-CT 


32 


2.5 


37.7 


236 


1210 


0.0 


Y 


a 


164n256vzw-rj 


64 


5.0 


302 


1890 


9660 


0.0 


Y 


cr 


164n256vzw-cr-nagb 


(51 


5.0 


302 


1890 


9660 


0.0 


N 


a 



a The 'vzw' suffix refers to the momentum-driven winds with Lr varying between 1.05-2.0, metallicity-dependent 
''wind! an d an extra kick to get out of the potential. 
b Box length of cubic volume, in comoving h _1 Mpc. 

c Equivalent Plummer gravitational softening length, in comoving h~ 1 kpc. 
d All masses quoted in units of 10 6 Mq. 
e Minimum resolved galaxy stellar mass. 



z = 1.5, despite containing a smaller volume unable to build 
larger structures. 



2.1 Metal Yields 

In previous Gadget-2 simulations including SH03b and 
OD06, metal enrichment was tracked with only one variable 
per SPH particle representing the sum of all metals and 
was assumed to arise from only one source, Type II SNe, 
which enriched instantaneously. While this is reasonably 
accurate when considering oxygen abundances, the abun- 
dances of other species can be significantly affected by al- 
ternate sources of metals. 

We have implemented a new yield model that tracks 
four species (carbon, oxygen, silicon, and iron) from three 
sources (Type II SNe, Type la SNe, and AGB stars) all 
with metallicity-dependent yields. These sources have quite 
different yields that depend significantly on metallicity, and 
inject their metals at different times accompanied by a large 
range in energy feedback. A more sophisticated yield model 
is required to model metal production from the earliest stars, 
abundance variations within and among galaxies, abundance 
gradients within the IGM, and abundances in the ICM. 

The four speci es chosen not only make up 78% of all 
metals in the sun (| Anders fc Grevessd Il989l ), but are the 
species most often observed in quasar absorption line spec- 
tra probing the IGM, X-ray spectra of the ICM, and the ISM 
of galaxies used to determine the galaxy mass-metallicity 
relationship. Furthermore, because these metals are among 
the most abundant, they are also often the most dynami- 
cally important when considering metal production in stars, 
the multi-phase ISM, and metal-line cooling of the IGM. 
We have not implemented metal-line cooling per individual 
species, but this may be straightforwardly incorporated in 
the future. 



2.1.1 Type II Supernovae 

Type II SNe enrichment follows that presented in SH03a, 
namely their equation 40 where gas particles are self- 
enriched instantaneously via 

AZspccics = (1 — /sN)j/spccics(Z):r— (1) 

where /sn is the fraction of the stellar initial mass func- 
tion (IMF) that goes supernova, x is the fraction of an 
SPH gas particle in the cold ISM phase, and t, is the 
star formation timescale. Our modification is that we fol- 
low the yield of each species individually using metallicity- 
dependent y ields, ?/ S pecies(^), from th e nucleosynthetic cal- 
culations by iLimongi fc Chiefnl |2005l ) instead of assuming 
y = Z Q = 0.02 as SH03b and OD06 did. Their grid of mod- 
els include SNe ranging from 13 to 35 Mq and metallici- 
ties ranging from Z = — 0.02. Using the total metallic- 
ity of a gas particle (i.e. the sum of the four species di- 
vided by 0.78 to account for other species), we employ a 
look up table indexed by meta llicity to obtain the 2/ sp ecies(Z). 
The lLimongi fc Chiefnl (|2005l ) yields are the most complete 
set of metallicity-dependent yields since IWooslev fc Weaverl 
1 19951 ). Both papers find similar yields for carbon and oxy- 
gen, the two spe cies most i mpor tant for IGM observations. 

We use the IChabriel (2003) IMF, although it is quite 
possible to modify the IMF in the future so as to have a 
top-heavy IMF under different conditions (e.g. |Pavell2008l ). 
We assume all stars between 10-100 Mq go supernova, com- 
prising /sn = 0.198 (i.e. 19.8% o f the total stellar mass in 
the IMF). We use the yields of ILimongi fc Chieffil (|2005T l 
comprising 13 — 35Mq over this larger mass range, thus as- 
suming the similar yields from stars between 10— 13Mq and 
35— IOOMq. Ot her supernova yield mo d els that include more 
massive stars (|Portinari et al.l Il998l ; iHirschi et all 120051 ) 
show higher carbon and oxygen yields from stars over 40 Mq 
at solar metallicity, but do not cover the range of metallici- 
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ties of lLimongi fc Chieffil (|2005h . The fraction of stars going 
supernova, /sn, is nearly twice as high since SH03a assumes 
a Salpeter IMF (/sn = 0.1), because there is a turnover 
at masses less than 1 Mq. The remaining 80.2% form star 
particles from which AGB feedback arises at later times. 



2.1.2 Type la Supernovae 

Type la SNe are believed to arise from mass accre- 
tion from a companion star that increases the mass 
beyond the Chandrasekhar limit, causing an explosion. 
Recently, the Type la SNe rate was measured by 
iMannucci et al.| (120051). and the resul t ing d ata was pa- 
rameterized by IScannapieco fc Bildstenl (120051 ) with a two- 
component model, where one component is proportional to 
the stellar mass ("slow" component), and the other to the 
star formation rate ("rapid" component): 

SWi?ia = AKL + BMi, M /yr (2) 

IScannapieco fc Bildstenl (l2005h determined that the best fit 
to the Mannucci et al.l (|2005l ^ data was provided by A = 
4.4 x 10" 14 yr" 1 and B = 2.6 x 10" 3 , with a time delay of 
0.7 Gyr in the slow component for the onset of Type la SNe 
production. 

To implement this in Gadget-2, we calculate the num- 
ber of Type la SNe formed at each timestep for every gas 
particle (from the first part of eqn. and every star parti- 
cle where the star was formed more than 0.7 Gyr ago (from 
the second term in eqn. [2} . Each Type la SN is assumed to 
add 10 51 ergs of energy, which is added directly to the gas 
particle or, in the case of star particles, added to the near- 
est gas particle. Each Type la is also assumed to produce 
O.O5M of carbon, 0.143M© of oxygen, 0.150M Q of silicon, 
and 0.613M Q of iron with these fou r species making up 78% 
of all metals made by Type la SNe l|Thielemann et ahll 19861 ; 
iTsuiimoto et al.lll995h . 



2.1.3 AGB Stars 

Feedback from AGB stars co mprise at least half of the to - 
tal mass returned to the ISM l|Wallerstein fc Knapp|ll99"^ ). 
AGB stars copiously produce carbon, referred to as car- 
bon stars, and other isoto pes of carbon, nitrogen, and oxy- 
gen (|Renzini fc Volilll98ll ) on a delayed timescale compared 
to the relatively instantaneous enrichment of Type II SNe. 
Heavier elements such as silicon and iron that remain un- 
processed by low-mass stars can now be returned back into 
the ISM instead of being trapped in stars. The mass and 
metallicity feedback from AGB stars is considerable and in 
some regions even dominant over supernovae. However, rel- 
ative to supernovae, the energy feedback can be considered 
negligible, since most of the mass leaves AGB stars at far 
less than 100 kms -1 . 

First we consider the mass feedba ck as a function of 
time for a star particle. We use the iBruzual fc Charlotl 
(2003) stellar synth esis models using the Padova 1994 
(jBressan et al.lll993h libraries of stellar evolution tracks to 
determine the mass loss rate from non-supernova stars as a 
function of age given a Chabrier IMF. Mass loss rates are 
calculated at an age resolution of 0.02 dex for six different 
metallicities (Z=0.0001, 0.0004, 0.004, 0.008, 0.02, and 0.05) 



covering stellar ages from log(t) = 7.18 — 10.14 (i.e. the age 
of the death of the most massive AGB stars to the age of 
the Universe). We interpolate in age and metallicity. 

To implement this in Gadget-2 we use each star par- 
ticle's age and metallicity to determine the mass loss rate 
from a lookup table generated from the lBruzual fc Charlotl 
|2003h population synthesis models. This is performed for 
every star particle each timestep with the total amount of 
mass lost being the product of the timestep and the mass 
loss rate. This mass is then transferred to the 3 nearest gas 
particle neighbors using a neighbor search. 

To illustrate the importance of mass feedback from low- 
mass stars consider a stellar population at Z — 0.02, which 
by log(t) = 7.34 yrs has returned about 17.5% of its mass via 
supernova feedback. By 100 Myr another 10.3% of the mass 
is returned to the ISM, and another 12.0% is returned by 1 
Gyr. More mass is returned to the ISM in the first Gyr from 
AGB stars than SNe. Another 9.7% is returned between 1 
and 10 Gyr. Overall, slightly greater than 50% of the mass 
of a Chabrier IMF is returned to the ISM with more than 
30% coming from low and intermediate mass stars. 

To model the retu rn of metals in to the ISM, we use the 
stellar yield models of iHerwid (I2004T) ( Z = 0.0001), iMarigol 
l|200ll ) (Z =0.008), and iGavilan et all (|2005l ) (Z =0.0126, 
0.0200, 0.0317) for a variety of stellar masses. We take the 
age at which a given stellar mass ends its life, from which we 
generate an interpolated lookup table of yields for a given 
species analogous to the mass loss lookup table. The incre- 
mental parcel of mass lost during each time step for each 
star particle is given a yield corresponding to the mass of a 
star dying at the star particle's age. This assumes that only 
one specific stellar mass is contributing to the entire mass 
loss from a stellar population at a given age, which is not 
a bad assumption considering that intermediate mass stars 
lose most of their mass during the short-lived AGB stage. 

We calculate AGB yield lookup tables for carbon and 
oxygen only. Silicon and iron remain almost completely un- 
processed in low and intermediate mass stars, so therefore 
we simply take their yields to be the original abundances of 
the star particle when it formed. To illustrate the yields, Fig- 
ure[T]plots the carbon yield as a function of age for a variety 
of metallicities. The dashed line is the relation between the 
Zero-age Main Sequence mass and the age of death. Carbon 
stars enrich copiously between ~200 Myr and 1 Gyr corre- 
sponding to stars of masses 2-4 Mq going through the AGB 
stage. The reason for this is that third dredge-up becomes 
efficient above 2 Mq transporting the products of double 
shell burning into the envelope (i.e. C 12 ) until hot-bottom 
burning becomes efficient above 4 -5 Mq transforming car- 
bon into nitrogen. iLia et al.l (|2002h also added delayed AGB 
feedback into SPH simulations, although they do not ex- 
tract yields directly from models; instead they use a time- 
dependent yield function without metallicity dependence for 
carbon. 

It is worth pointing out that our simulation do not track 
the metals in stellar remnants. The metal products arising 
from nucleosynthesis in AGB stars for the most part remain 
in the white dwarf remnant once it has blown its envelope off, 
and we do not track the creation of these metals in the star 
particles. The same is true for neutron stars and stellar black 
holes. Fortunately these metals remain locked in their stellar 
remnants on timescales far beyond the age of the Universe, 
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Figure 1. Carbon yields of gas lost from AGB stars for sev- 
eral metallicit i es calc ulated bv lHerwigl ||2004P ) . iMarigo l|200lh and 
iGavilan et al. I d2005l> converted from a function of stellar mass 
to a function of age by using a star's lifetime (dashed line). In 
Gadget-2, we take a star particle's age and assign the yield from 
these curves to the fractional amount of mass lost from the star 
particle and added to the 3 nearest SPH neighbors. The carbon 
yields jump between ~ 200 Myr and 1 Gyr when AGB stars in 
the range of 2-4 Mq are dying as carbon stars; AGB yields are 
highly dependent on age as well as metallicity. The Z = 0.008 
yields show the strongest yields, especially toward lower mass, 
while by Z = 0.0317 the strongest yields come from 5Mq stars. 



hence we can ignore them as a component of observable 
metals. When we talk about the global metal budgets in Sj4l 
we do not include metals trapped in stellar remnants. 



2.2 Group Finder 

A group finder added to Gadget-2 allows us to add new 
dynamics based on the properties of a gas or star particle's 
parent galaxy. This is especially important for momentum- 
driven winds, which MQT05 argued depend on the proper- 
ties of a galaxy as a whole (specifically a galaxy's a). We 
will show in H2.3I that winds derived from the group finder 
more accurately determine wind speeds for the momentum- 
driven wind model versus using the potential well depth of 
a particle as a proxy of galaxy mass as was done in OD06 
and DO07. 

Our group finder is based on the friends-of-friends 
(FOF) group finder kindly provided to us by V. Springel, 
which we modified and parallelized to run in situ with 
Gadget-2. Gas and star particles within a specified search 
radius are grouped together and cataloged if they are above 
a certain mass limit, which we set to be 16 gas parti- 
cle masses. The search radius is set to 0.04 of the mean 
interparticle spacing corresponding to an overdensity of 
(1/0.04) 3 = 15625. We chose this value after numerous com- 
parisons with the more sophisticated (but far more compu- 



tationally expensive) Spline Kernel Interpolative DENMAX 
(SKIDJj group finder performed on snapshot outputs. 

The FOF group finder is run every third domain decom- 
position, which corresponds to an interval usually between 2 
and 10 Myr depending on the dynamical time-stepping for 
the given simulation (i.e. smaller timesteps for higher res- 
olution). This timestep was chosen to be smaller than the 
dynamical timescales of most if not all galaxies resolved in a 
particular simulation. The galaxy properties output include 
gas and stellar mass and metallicities, star formation rate, 
and fraction of gas undergoing star formation. The addi- 
tional CPU cost ranges from 7-10% for a 64-processor run 
to less than 5% for a 16-processor run. Once the FOF group 
finder is run, each particle in a galaxy is assigned an ID so 
that it can be easily linked to all the properties of its parent 
galaxy. 

A comparison plot of the two group finders in Figure [2] 
for the 132n256vzw-cr simulation at z = shows very good 
but not perfect agreement for total baryonic galaxy mass 
(Afgalj Panel (a)). Galaxies are matched up by requiring a 
< 5fr, _1 kpc difference between SKID and FOF positions. 
Stellar masses (M», lower left panel) are nearly identical, 
but the associated gas mass (M gas , upper right) shows more 
scatter. Note that there is no explicit density or temper- 
ature threshold for including gas in the FOF case, but in 
the SKID case only gas with overdensities > 1000 compared 
to the cosmic mean are included; this may contribute to 
some of the scatter. FOF group finders have the tendency 
to group too many things together in dense environments, 
and this is most noticeable in the associated gas masses. 
There are significantly better agreements in gas masses at 
higher redshift, where dense group/cluster environments are 
less common and gas fractions are greater. 

The mass functions of both group finders when all 
galaxies are included show very good agreement (lower right 
panel of Figure [2}. SKID and FOF each find 6 galaxies with 



M ea i > 10^M©, 112 and 113 galaxies with 10^ > M ga i > 



10 xl M(s, 793 and 719 galaxies with 10 11 > M ga i > 1O 1O M , 
and 4296 and 3596 10 10 > Af ga i > 1O 9 M respectively. 
1.2 x 10 9 Mq is the galaxy ma ss resolution limit d efined as 
the mass of 64 SPH particles i|Finlator et al.ll2QQ6l ). Owing 
for the tendency of FOF to over-group satellite galaxies with 
central galaxies in dense environments, there is a deficiency 
of small satellite galaxies in the FOF case. The total amount 
of mass in all resolved galaxies is within 0.2% between the 
two group finders, however there is 24% more mass grouped 
to the 6 largest galaxies in the FOF finder. 

As a side note, our group finder has the flexibility to en- 
able modeling of merger-driven or mass-threshold processes 
such as AGN feedback. Although our simulat ions implic- 
itly i n clude both hot and cold-m ode accretion l|Keres et al.l 
l2005l) . |Pekel fc Birnboiml (|2006l ) suggests AGN feedback af- 
fects only hot mode accretion, and hence there exists a 
threshold halo mass above which AGN feedback is effec- 
tive. Our group finder can identify galaxies where AGN feed- 
back may be necessary to curtail star formation and drive 
AGN winds. Conv ersely, if AGN feedback is driven by th e 
onset of a merger (|Di Matteo. Springel fc Hernquistll2005l ). 
our group finder can identify mergers and add merger-driven 
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Figure 2. The comparison of two group finders run on a 32 h~ x Mpc 2 X 256 3 simulation at z = 0. FOF is the much faster algorithm 
run in situ with Gadget-2 at frequent intervals to determine galaxy properties on the fly, while SKID is a post-processing DENMAX 
algorithm run on simulation snapshots. Panels (a), (b), and (c) compare the masses (total, gas, and stellar mass respectively) of groups 
found within 5 kpc of each other by the two group finders. Stellar masses agree very well while the gas masses (a smaller component by 
z = 0) show more scatter. The total mass functions (Panel (d)) of both group finders lie nearly on top of each other above the resolution 
limit (dashed line), although FOF tends to group some satellites found by SKID with the largest few galaxies. The M ga i functions do 
not resemble a Schechter function, although when one plots M» there is a more pronounced turnover below M* = 10 10 AfQ. 



AGN feedback that curtails star formation. We leave such 
implementations of AGN feedback for future work, though 
we note that the heating of gas may transport a non-trivial 
amount of energy and metals into the IGM. 

Frequent outputs of the group finder allow us to track 
the formation history of galaxies to see how galaxies build 
their mass (e.g. through accretion or mergers). We can trace 
the star formation rate during individual mergers to see if 
our simulations are producing bursts of star formation. Out- 
putting group finder statistics during a run allows us to im- 
mediately look at integrated properties such as the galaxy 
mass function, the mass-metallicity relationship, and the 
specific star formation rate as a function of mass just to 
name a few. This is valuable as these relations may exhibit 



changes on timescales shorter than the simulation snapshot 
output frequency. 



2.3 Wind Model Modifications 

We use the same kinetic wind implementation of SH03a, 
whereby particles enter the wind at a probability rj, the mass 
loading factor relative to the star formation rate. Wind par- 
ticles are given a velocity kick, w w i n d) in a direction given by 
v x a (ie. perpendicular to the disk in a disk galaxy) . These 
particles are not allowed to interact hydrodynamically until 
either they reach a SPH density less than 10% of the star 
formation density threshold or the time it takes to travel 
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30 kpc at u w i n< j; the first case overwhelms the instances of 
the second case in our simulations. We admit that this phe- 
nomenological wind implementation insufficiently accounts 
for the true physics of drivi ng superwinds as wel l as the 
multi- phase aspect of winds ^Strickland et al.ll2002l ; iMartinl 
l2005bl i: see iDalla Vecchia fc Schavel (|2008l ) for an in-depth 
discussion of some of the insufficiencies of such winds in sim- 
ulations. However we want to stress that while we cannot 
hope to model the complexities of the outflows, the focus of 
this paper primarily depends on the much longer period of 
subsequent evolution. We use the momentum-driven wind 
model with variable w w i n d and rj, due to its successes in pre- 
vious publications mentioned earlier. 

In the momentum-driven wind model analytically de- 
rived by MQT05, w w i n d scales linearly with the galaxy ve- 
locity dispersion, ct, and r\ scales inversely linearly with ct. 
We again use the following relations: 



"wind 

V 



£0 

3 

a 



(3) 
(4) 



where ft is the luminosity factor in units of the galactic 
Eddington luminosity (i.e. the critical luminosity necessary 
to expel gas from the galaxy), and cto is the normaliza- 
tion of the mass loading factor. The outflow models used 
in this paper are all of the 'vzw ' variety described in §2.4 of 
OD06, in which we randomly select luminosity factors be- 
tween fx, — 1.05 — 2.00, include a metallicity dependence for 
/l owing to more UV photons output by lower- metallicity 
stellar populations 



h = /l,© x 10 



-0.0029*(log Z+9) 2 - 5 +0. 417694 



(5) 



and add an extra kick to get out of the potential of the 
galaxy in order to simulate continuous pumping of gas until 
it is well into the galactic halo (as argued by MQT05). 

One difference is that we use cto = 150 kms -1 instead 
of 300 kms -1 in the equation for r\ — uo/cxSFR. Since our 
assumed WMAP-3 cosmology produces less early structure 
that the WMAP-1 cosmology used in OD06, it requires less 
suppression of early star formation and hence lower mass 
loading factors (DO07). We find that <xo = 150 kms -1 with 
the WMAP-3 cosmology generally reproduces the successes 
of cto = 300 kms" 1 with the WMAP-1 cosmology. 

The main modification we make to the outflow imple- 
mentation is how ct is derived. Putting the wind parameters 
in terms of ct is the most natural, because the fundamental 
q uanti ty MQT05 use to derive momentum-driven winds is 
\J £M, which equals \/2ct for an isothermal sphere. In pre- 
vious runs without a group finder , we derived ct using the 
virial theorem where ct = \J—\<&, with $ being the gravi- 
tational potential at the initial launch position of the wind 
particle. We will call this old form of the wind model <E>- 
derived winds. While this derivation of ct should adequately 
work for an isothermal sphere, the <E> calculated by Gadget- 
2 is the entire potential: the galaxy potential on top of any 
group/cluster potential. As galaxies live in groups and clus- 
ters more often low redshift, the wind speeds from "^-derived 
winds become overestimated. To counteract this trend, we 
artificially implemented a limit of ct = 266 kms -1 , which 
corresponds to a M ga i ~ 2.7 x 10 12 M@, the mass of a giant 
elliptical galaxy at z = 0. In the deep potential of a clus- 
ter, all galaxies no matter what size would drive winds at 



the speed of this upper limit. The overestimate of "wind pre- 
vented us from trusting our wind model at lower redshifts; 
therefore we usually stopped our simulations in OD06 at 
2 = 1.5. 

We now introduce er-derived winds, where ct is calcu- 
lated from Mgai as determined by the group finder. We 
use the same r elation as MQT05 (their equation 6) from 
iMo et alJ 1 19981 ) assuming the virial theorem for an isother- 
mal sphere: 



ct = 200 



M ga l Sl m H{z) 



5 x 10 12 Q. b 



Ho 



1/3 



kms 



(6) 



Since our group finder links only baryonic mass, we multiply 
Mgai by flm/^lb to convert to a dynamical mass. The wind 
properties are hence estimated from the galaxy alone. We do 
not calculate ct from the velocity dispersion of star and/or 
gas particles in the galaxy, since our tests show the resolution 
is insufficient to derive a meaningful ct. 

The H(z)/Ho factor in equation [6] increases ct for a 
given mass as H(z) increases at higher redshift. For exam- 
ple, H(z = 6) = 10.2 x Ho resulting in winds being 2.17x 
as fast being driven from the same mass galaxy at z = 6 
versus z = 0. Physically, galaxies that form out of density 
perturbations at higher z have overcome faster Hubble ex- 
pansion, and therefore their ct is higher. The same mass 
galaxy formed at higher redshift with a higher ct means 
the higher redshift galaxy must be more co mpact. Such a 
scenar io is supported by t he ob servations bv lTruiillo et alJ 
l|2006l ) and iTruiillo et all l|2007h showing a trend of galax- 
ies b ecoming smaller from z — — > 2, in general agreement 
with lMo et ail (|l998l ). A more compact galaxy drives a faster 
wind in the momentum-driven wind model, because M«md is 
a function of M ga i/r, not just M ga i; there is more UV pho- 
ton flux impinging each dust particle. We will show that 
increasing "wind emanating from the same mass galaxy to- 
ward higher z has important consequences in enriching the 
IGM at high-z while not overheating it at low-2. 

Another modification to the wind model we add are new 
"wind speed limits. The first depends on the star formation 
timescale, tsfr. The momentum-driven wind equation de- 
rived by MQT05 and used in OD06 and DO07 assumes a 
starburst occurs on the order of a dynamical timescale, td, 
which is often the case in a merger-driven starburst. How- 
ever, MQT05 also derive a maximum ct, CT ma x, above which 
a starburst cannot achieve the luminosity needed to expel 
the gas in an optically thick case. Although it is not clear 
how CT max varies with the tsfr, we modify their equation 23 
to have an inverse linear dependence on the star formation 
timescale 



4000 x kms" 1 
tsfr 



(7) 



and assume a td of 50 Myr. The end result is a reduction 
of 5-10% in the average " W i n d out of M ga i ~ 10 12 at z < 1. 

A second speed limit we impose allows no more than 2x 
the total SN energy to be deposited into the wind. This does 
not violate the energetics as the energy for these winds is 
coming from momentum deposition over the entir e lifetime 
of a s tar, which is of the order 100 x the SN energy l)Schaererl 
2003). This limit was instituted to disallow extreme values of 
"wind (i.e. > 1500 kms -1 ) emanating from the most massive 
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l°g(M gal ) 

Figure 3. Average launch wind speed as a function of M ga i (de- 
rived using SKID) versus dotted lines for the predicted (t> w i n d) 
assuming no speed limits. The wind speeds, which are derived us- 
ing the FOF finder during a run, agree very well with the masses 
determined by SKID in the post-processing stage. The decline in 
(■"wind) with redshift for a given M ga i is a result of the H (jz) de- 
pendence in equation [6] This trend loosely matches observations 
of faster winds for a given mass galaxy at high redshift. 

haloes, reducing u w ind at most 30-50% in the most massive 
galaxies at z < 1. 

In Figure [3] we plot the average wind speed, (v W md) , 
as a function of galaxy mass at four redshifts for a variety 
of box sizes ranging from 8 to 64 fo -1 Mpc along. Dotted 
lines show the predicted (« w ind) for solar metallicity, assum- 
ing no speed limits. As with every plot in this paper, M ga i 
is the SKID-derived baryonic mass, not the FOF-derived 
mass from which Wwmd is calculated. The simulations repro- 

1/3 

duce the predicted trend t> W md oc M ^ . Divergences at low 
mass result from faster winds being driven by low-mass, low- 
metallicity galaxies as well as some satellites being grouped 
with a central galaxy in groups/clusters by FOF. The lat- 
ter effect appears to be sub-dominant though as evidenced 
by the overlap of the relation among simulations at differ- 
ent redshifts. The deviations at the high-mass end at all 
redshifts are dominated by the second wind speed limit dis- 
cussed above. 

The reduction in v w ind for a given mass galaxy due to 
the H(z) dependence in equation [5] is the reverse of the 
trend in the "^-derived models, and are in better agreement 
with low-2 observations. In the new cr-derived wind model, 
a lO n M0 galaxy (such as the Milky Way) launches an av- 
erage wind particle at 790 kms -1 wind at z — 3 and 450 
kms -1 at z — while the corresponding values for the $- 
derived winds without any speed limits are 1040 kms -1 and 
1220 kms -1 respectively. The la tter are far abo ve the val- 
ues observed in the local Universe. iMartir] (|2005al ) found the 
relation Utorm oc 2.1 x v c i IC fit her observations best, where 
utcrm is the terminal velocity of the wind. For a 10 n Mo 
galaxy, this corresponds to v c i IC ~ 125 kms -1 leading to 
utcrm ~ 265 kms -1 . This is much nearer our cr-derived 
value once the extra velocity boost to leave the potential 
of the galaxy is subtracted. High-velocity clouds (HVC's) 



may be material that is blown into t he halo from the Milky 
Way l|Wakker fc van Woerden|[l997l ). and these have veloc- 
ities that are <C 1000 kms -1 , in better agreement with 
velocities expected in the cr-derived model. Of course the 
Milky Way is not a star-bursting galaxy driving a power- 
ful wind; however, it is still possible that it is kicking up a 
significant amount of gas into the halo. Indeed, as we will 
discuss later ( £|5.4[) . the outflows in our simulations don't al- 
ways reach the IGM, but particularly at low-2 may be more 
aptly described as "halo fountains" , where the outflows only 
propagate out to distances comparable to the galactic halo. 

It remains very difficult to determine from observations 
whether wind speeds increase for a given mass galaxy at 
high-z as the cr-deri ved winds pred i ct. O bserving asymmet- 
ric Lyman-a profiles fWilman et al.l (|2005T ) find what they in- 
terpret as 290 kms -1 outflows arou nd a LBG with 1 11 M ( y ) 
baryons at z = 3.09. More recently ISwinbank et alj (|2007T ) 
observe outflows up to 500 kms -1 at z = 4.88 around a 
lensed galaxy with a dynamical mass as low as 10 10 M Q . 
These observed outflows appear to be beyond the virial ra- 
dius in both cases and correspond to our velocities once we 
have subtracted the extra velocity bo ost we add to get out 
of the potential well. Our wterm for the lWilman et alT(|2005l ) 
z = 3 galaxy would be closer to 500 kms -1 , somewhat 
above their observed values. We expect f w j n d ~ 300 kms -1 
(wterm ~ 200 kms -1 ) for the lSwinbank et al.l (|2007h object, 
but their mass is only a lower limit. If this object is a 1O 1O M0 
baryonic-mass galaxy, we would derive wtcrm ~ 400 kms -1 . 
Overall, the cr-derived winds predict velocities that are at 
least in the ballpark of observed values. Additionally, sur- 
veys of LBG's atz ~ 3 l|Pettini et al.ll200ll ; IShaplev et~all 
120031 ) find outflows of several hundred kms -1 to be ubiq- 
uitous, often driving an amount of mass comparable to the 
star forming mas (i.e. 77 ~ 1) 



3 THE UNIVERSAL ENERGY BALANCE 

Armed with these new simulations, we can now investigate 
how outflows move mass, metals, and energy around the Uni- 
verse. In this section we focus on the energetics of outflows, 
and its impact on cosmic star formation and temperature. 
We will compare the new cr-derived wind model's behavior 
versus the old "I>-derived wind model, and study the impact 
of AGB feedback. Specifically, cr-derived winds inject much 
less energy at late times making a cooler and less-enriched 
IGM while leading to more star formation. The inclusion 
of AGB feedback does not really affect the global energy 
balance, but does increase the number of stars formed and 
more significantly affects the amount and location of metals, 
as we discuss in 

The his tory of the cosmic star formation rate den- 
sity (SFRD; iMadau et aT]|l996l ; IlTiIv et all 119961 ) is a key 
observabl e that has received much attention in r ecent 
years (e.g. lHopkins fc Beacom 2006; Fa rdal et al.ll2007l ). The 
SFRD plot in the upper l eft panel of Figure [4] show s how 
our models compare to the lHopkins fc B eacom (2006) com- 
pilation (black lines, two different fits). It should be noted 
that we are resolving galaxies down to « 1O 1O M0 in the 
models shown with solid lines, and hence the star forma- 
tion density at z > 3 is increasingly underestimated (see 
e.g. SH03b). Hence this should be considered as an illustra- 
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Figure 4. Three 32 ft - 1 Mpc 2 X 128 3 simulations: <E>-derived winds including AGB feedback (green), cr-derived winds with AGB 
feedback (magenta), and cr-derived winds without AGB feedback (blue) demonstrate how global quantities evolve as a result of energetic 
interactions between star formation-driven feedback from galaxies and the state of the IGM. Panel (a) shows the star formation density 
histories compared to the | H o D kins fc Beaco ^ B compilation averages (black lines for segmented and non-segmented fit, scaled to a 
Chabrier IMF). Panel (b) shows the average virial ratio with which wind particles are launched. Solid/dashed lines in Panel (c) trace the 
volume/mass- weighted average gas oxygen content, and Panel (d) shows the volume- weighted temperature. The new cr-derived winds 
do not inject as much energy or metals into the IGM while increasing the star formation. AGB feedback has a relatively minor effect 
in the energetic balance between galaxies and the IGM. The lower resolution 132nl28 simulation misses about half of all star formation, 
therefore we also show the 132n256 (dash-dotted line) and the 116n256 (dotted line) simulations in Panel (a), which both agree better 
with observations. 



tive comparison, whose main point is to highlight the differ- 
ences between various models. A more complete SF history 
at early times is shown by the dashed and dotted magenta 
lines from our higher resolution 2 x 256 3 particle simula- 
tions with cr-derived winds with AGB feedback. The metal 
enrichment of the IGM at early times is significantly higher 
when the star formation from smaller galaxies is included; 
therefore we use only higher resolution models to explore 
the IGM enrichment at z > 2. Also, our models suggest an 
overestimated SFRD at the lowest redshifts. This indicates 
that some other form of feedback is needed to suppress star 



formation in massive galaxies at l ate times (such as AGN 
feedback; e.g. ICattaneo et al.ll2007r ). 

The three models shown, the ^-derived winds and the 
cr-derived winds with and without AGB feedback, are indis- 
tinguishable in their global SFRD's above z — 4 because the 
mass loading factor and not wind speed is the largest deter- 
minant of star formation (OD06). As explained in OD06, the 
earliest wind particles have not had time to be re-accreted 
by galaxies, therefore the SFR is regulated purely by how 
much mass is ejected (77). The faster wind velocities of the 
^-derived model at z < 2 decrease star formation relative 
to the cr-derived model, because wind particles are sent fur- 
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ther away from galaxies making this gas harder to re-accrete 
while also heating the IGM more and further curtailing star 
formation. We will quantify this recycling of winds in 

The average virial ratio of winds (-E ra t), defined as the 
ratio of the kinetic energy to the potential at the launch posi- 
tion (0.5 x u w i n d 2 / — <£"), in the upper right panel shows how 
cr-derived winds inject progressively less energy into their 
surroundings with time. E ra x should be invariant across time 
for ^-derived winds in its simplest form, but falls sharply be- 
low 2 = 1 due to our wind speed limits, and rises slightly 
at high-z due to the metal dependent « W md- The faster 
derived winds spatially distribute more metals (solid lines 
in lower left panel) and heat the IGM (lower right panel) 
to a greater extent than the cr-derived winds. We will show 
in ij4.2l that the cooler, less-enriched IGM of the cr-derived 
wind models makes a significant difference in metals seen 
in quasar absorption line spectra. The gaseous metallicity 
(dashed lines in lower left panel) is slightly higher in the 
^-derived wind model despite fewer overall stars formed, 
because fewer metals end up in stars. 

Finally, we consider the impact of AGB feedback. AGB 
stars do not add any appreciable energy feedback, as demon- 
strated by the invariance in the _E rat and volume-averaged 
temperature in the models with (magenta lines) and without 
(blue) AGB feedback. However, AGB stars provide feedback 
in the form of returning gas mass to the ISM, which is now 
available to create further generations of stars. For instance, 
with AGB feedback the SFRD at z — is increased by nearly 
a factor of two compared to without. Hence models that do 
not include such stellar evolutionary processes may not be 
correctly predicting the SFRD history. 



4 THE UNIVERSAL METAL BUDGET 

In this section, we will investigate what the sources of cos- 
mic metals are, and where these metals end up. To do so, we 
examine quantities summed over our entire simulation vol- 
umes, with special attention to the three different sources of 
metals: Type II SNe, Type la SNe, and AGB stars. We also 
discuss the impact of AGB feedback and cr-derived winds. 



4.1 Sources of Metals 

4-1.1 The Stellar Baryonic Content 

Since stars produce metals, we first examine the evolution of 
the stellar baryonic content. The key aspect for metal evolu- 
tion is that stellar recycling provides new fuel for star forma- 
tion. For the Chabrier IMF assumed in our simulations, su- 
pernovae return ~ 20% of stellar mass instantaneously back 
into the ISM (which is double that for a Salpeter IMF), and 
delayed feedback will eventually return another ~ 30% over 
a Hubble timescale. With half the gas being returned to the 
ISM, most of it (~ 80%) on timescales of less than a Gyr, 
subsequent generations of stars can form, leading to greater 
metal enrichment. The difference between a-derived wind 
models with and without AGB feedback in the SFRD plot 
(upper left panel of Figure [4| demonstrates how AGB feed- 
back makes available more gas for star formation at later 
times. 

Despite more star formation in the 164n256vzw-cr model 
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Figure 5. Global quantities binned in 200 Myr bin of stars formed 
(<5f2«, solid black line) and mass lost via delayed feedback (<5f2delj 
dashed red line) relative to Q b in the 32 fc _1 Mpc 2 X 256 3 simu- 
lation. The ratio of the two is the long-dashed blue line. 



with AGB feedback, slightly more mass is in stars at z — 
in the no- AGB feedback model by z = (Q*(z = 0)/Qb = 
0.081 vs. 0.071) since there is no mass loss from long-lived 
stars. However, if we count all stars (including short-lived 
stars undergoing SNe) formed over the lifetimes of our sim- 
ulated universes, the AGB feedback model forms 37% more 
stars (tl,(all)/Slb = 0.138 vs. 0.101). This is the more rele- 
vant quantity when considering the metal budget of the local 
Universe. Another way to think of this is that the stars in to- 
day's local Universe account for (£l*(z = 0)/fi*(aZZ)=) 52% 
of all the stars that have ever existed, assuming the mass 
in short-lived stars is negligible (a safe assumption in the 
low-activity local Universe). 

Figure[S]shows the amount of baryons formed into stars, 
<50» (i.e. the star formation rate), and the amount of stellar 
mass lost via delayed feedback, 5fidelj as a function of time. 
The ratio SQ^d/S^* increases as more generations of stars 
are continuously formed with each generation contributing 
to 5f2del' The quantity 8£ld c i/SQ* should approach 0.3 for a 
steady star formation rate in a Hubble time, because 30% 
of the stellar mass is returned via delayed feedback, but this 
is actually exceeded since the star formation is declining at 
late times. The amount of material lost via delayed feedback 
over the history of the Universe is found by subtracting the 
stars at z — from the number of long-lived stars formed 
over the history of the Universe. 
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14Gyr 



tffi*(t)(l - /sn)«K - n*(2 = 0) 



(8) 



OGyr 



where we assume fi* (z — 0), the stars in today's Universe, 
has a negligible quantity of short-lived stars. The amount 
of delayed feedback is fidei = 0.039 or 55% of the stellar 
baryons in today's Universe, which is a significant quantity. 
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4.1.2 The Metal Content 

We segue into our discussion of metals by plotting the global 
production of the 4 species tracked by their source (Type 
II SNe, Type la SNe, or AGB stars) in the left panels of 
Figure Metals produced by Type II SNe depend on the 
current star formation rate (<5f2*), which is apparent by the 
fact that the solid lines in the left panels have roughly the 
similar shape to the global stellar mass accumulation rate 
(Figure[5|. Type II SNe dominate the enrichment for all four 
elements at all redshifts, hence global chemical enrichment 
is reasonably approximated by current star formation alone. 

Dividing the amount of each type of metal formed via 
Type II SNe by f2„ gives the SN yield of that element shown 
as solid lines in the right panels of Figure |S] These yields 
should and do match the Type II SNe yield tables that are 
an input to the simulations. The yields do not evolve sig- 
nifica ntly because there is only weak metallicity dependence 
in the lLimongi fc Ch icffi (2005j) yields employed here. Note 
that previous work has generally assumed no metallicity de- 
pendence. Above z = 6, the yields are slightly lower due to 
less metals injected by low-metallicity stars. A slight upturn 
is noticeable at lower redshift for carbon and silicon as their 
yields increase with metallicity. 

Turning to metals injected via AGB feedback (dashed 
red lines in Figure |BJ), we find more complex behavior that 
varies among the species. The left panels show that AGB 
feedback is an important source of carbon (~ 30% by z — 0), 
iron and silicon (both ~ 25% by z = 0), but is a negligible 
contributor for oxygen. Dividing these values by fidei results 
in the global AGB yields as a function of redshift in the right 
panels. For carbon and oxygen, these are a summation of the 
yield values in our input tables for all stars of various ages 
and metallicities undergoing AGB feedback. 

The carbon yields shows the most interesting evolution, 
which are indicative of the underlying processes of stellar 
evolution in different stars. As explained in H2.1.3l and shown 
in Figure[T] the most massive and short-lived AGB stars (4-8 
Mq) have hot-bottom burning that destroys carbon; these 
are the stars losing mass via AGB feedback the most at 
very high-z. Between 2-4 Mq, the third dredge-up makes 
AGB stars into carbon stars with very high carbon yields 
from stars dying 200 Myr to 1 Gyr after their formation. At 
z ~ 2, carbon stars dominate the carbon yields, but then less 
massive stars (< 2M©) without the third dredge- up begin to 
reach the AGB phase, and the ensemble AGB carbon yield 
begins to decline. 

Oxygen is burned in AGB stars, resulting in a net de- 
crease in its overall content as a result of delayed feedback. 
Accounting for a minor contribution from Type la SNe, the 
vast majority of gaseous oxygen is synthesized in Type II 
SNe. Hence oxygen is the ideal species to trace the cosmic 
evolution of Type II SNe. 

The AGB yields of silicon and iron may at first be sur- 
prising considering that AGB stars do not process these ele- 
ments. These yields reflect the ensemble metallicities of mass 
loss from AGB stars, since they neither create nor destroy 
heavier elements at any significant rate, but instead sim- 
ply regurgitate them. Most surprising is that more iron is 
ejected from AGB stars than Type la SNe. The difference 
between these two forms of feedback associated with stars 
is that the iron yield of Type la SNe is nearly a half (i.e. 



0.6 Mq formed per 1.4 Mq SNe) and should significantly 
enrich its local environment, while the iron lost from AGB 
stars should have a slightly lower yield than the surround- 
ing gas metallicity since these stars are older and hence less 
enriched. 

It is curious and probably not correct that iron and 
silicon AGB yields exceed solar metallicities by as much a 
factor of 2 — 3 by z = 0. This means that at z = 0, the 
average AGB star is at least 2 — 3 x Zq, which is almost 
definitely too high when s tars younger than the Sun in the 
Milky Way disk are ~ Zq l |Twaroelll98"oh . Even though most 
z — AGB mass loss comes from stars younger than the Sun 
since most AGB feedback occurs within 1 Gyr for a Chabrier 
IMF, the extremely super-solar metallicities are indicative of 
too much late star formation. Reasons for this include too 
much star formation in massive systems in our simulations, 
as well as our slightly high value of Qb = 0.048 (instead of 
the currently more canonical Qt ~ 0.044). 

Table [2] summarizes the sources of metals at z = 2 
(roughly 10 10 years ago) and z = for the 132n256vzw- 
cr simulation. These can be compared to available observa- 
tional constraints. Using oxygen, the global metallicity av- 
eraged over all baryons is (Zb)(z — 2) ~ 0.064Z© rising to 
{Zh)(z — 0) ~ Q.23/?0 . These value s are r emarkably simi- 
lar to those derived bv lBouche et al.l (|2007T ) (hereafter B07) 
({Z b )(z = 2) ~ 0.056 and (Z b )(z = 0) ~ 0.20Z©) where 
they assumed a Salpeter IMF-weighted metallicity yield of 
0.024 and inte grated over t he sta r formation history of the 
Universe from Icole et~ai1 (|200ll ). While encouraging, this 
comparison is highly preliminary owing to many systemat- 
ics, such as the fact that our simulations produce too many 
stars overall, and the assumption of a Salpeter IMF at all 
times is pr obably not consistent with observations (see e.g. 
|Pavel200a . and discussions therein). It is hoped that improv- 
ing observations will enable more interesting constraints on 
cosmic chemical evolution models. 

While it is well-known that Type II SNe dominate car- 
bon, oxygen, and silicon production, it may be somewhat 
surprising for iron, considering that Type la SNe are often 
assumed to be the primary producers of iron; however, this 
is actually only true in environments dominated by older 
stars. Long-lived stars also destroy oxygen, eliminating 20% 
of the oxygen formed by Type II SNe. Processing of oxy- 
gen by AGB stars helps to move oxygen abundances from 
alpha-enhanced levels to solar levels. Of course, long-lived 
stars do make a net surplus of both carbon and oxygen in 
post-Main Sequence nucleosynthesis, however most of these 
metals remain locked in stellar remnants, which we do not 
track in our simulat i ons a nd are not included in this table. 
iFukugita fc Peebles! l|2004l ) estimate a (Z b )(z = 0) ~ O.68Z 
for all metals including those in remnants, which exceeds the 
metals not locked up (i.e. the ones we track) by a factor of 
a few. 

Even though AGB feedback injects an appreciable 
amount of carbon into surrounding gas (~ 40% of carbon 
injected via Type II SNe), the net surplus of carbon result- 
ing from AGB feedback is only ~ 5% of Type II SNe by 
z — 0, because much of this carbon is coming from stars 
with high metallicity already. Carbon stars (~ 2 — 4Mq) 
add to the overall cosmic carbon abundance while higher 
and lower mass AGB stars reduce the amount of carbon. A 
larger impact on carbon production comes from recycled gas 
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Figure 6. Metal production by species relative to f2b along with their yields plotted in 200 Myr bins in the 32 h~ 1 Mpc 2 x 256 3 
simulation. Panels (a-d) trace the amount of carbon, oxygen, silicon, and iron respectively returned to the gas phase by Type II SNc 
(black solid lines), Type la SNe (green dotted lines), and AGB stars (red dashed lines). The yields for the various species are calculated 
in Panels (e-h) by taking lines in Panels (a-d) and dividing by the matching line types in Figure [5] long dashed lines on the left indicate 
lAnders fc Grcvcsse (1989) solar values. Type II SNe yields remain relatively constant despite metallicity-dependent yields. The AGB 
carbon yields strongly depend on z and peak when 2-4 Mq die. AGB silicon and iron yields are simply reprocessed metals formed in 
SNe and reflect the ensemble metallicities of all stars undergoing AGB feedback at a given redshift. 
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Table 2. Sources of Feedback 



Species 


n/n b 


% Type II 


% Type la 


% AGB Stars 


% AGB Processing 


2 = 2 


Baryons 




0.87 


.0016 


0.83 




Carbon 


1.677e-04 


85.70 


0.0" 


23.25 


+14.30 


Oxygen 


6.123e-04 


105.4 


0.32 


3.48 


-5.75 


Silicon 


5.642e-05 


96.50 


3.57 


9.30 


+0.00 


Iron 


5.013e-05 


83.49 


16.51 


9.23 


+0.00 


2 = 


Baryons 




3.74 


.0098 


5.42 




Carbon 


6.566e-04 


95.35 


0.0° 


38.29 


+4.65 


Oxygen 


2.194e-03 


126.1 


0.52 


7.68 


-26.60 


Silicon 


2.549e-04 


95.29 


4.72 


33.89 


+0.00 


Iron 


2.265e-04 


78.32 


21.68 


33.16 


+0.00 



"These models were run with no Type la yields for carbon. These carbon yields are fractionally insignificant anyway. 



that enables more Type II SNe; as noted before, 37% more 
stars form in simulations that include AGB feedback. Stars 
also lose mass via AGB feedback when they have moved 
away from the sites of their formation, and can directly en- 
rich metal-poorer areas such as the ICM and intra-group 
medium. Hence the location of feedback from AGB stars 
and Type la SNe turns out to be important for understand- 
ing enrichment in various environments. We consider this 
topic next. 



4.2 The Location of Metals 

B07 calculated from observations that metals migrated from 
gas to stars between z = 2 — > as metals fall back into 
the deeper potential wells of growing galaxies, and are more 
likely to remain there as star formation-driven winds decline 
at low-z. Our new simulations generally agree with the re- 
sults of B07 that about one-third of metals are in stars at 
z = 2, increasing to two-thirds by z = 0. 

In Figure [7] we subdivide mass and metals further by 
their baryonic phases. The top panel shows the baryonic 
fraction in diffuse gas (T < 10 ), warm-hot intergalactic 
medium (WHIM) gas (10 5 < T < 10 7 ), hot gas in clus- 
ters (T > 10 7 ), and galaxies (stars and ISM) for both the 
164n256vzw-cr, 164n256vzw-a-nagb, 132nl28vzw-<I> models as 
a function of redshift. The steady transfer of gas from dif- 
fuse to WHIM and hot IGM b etwe en z = 3 -> wa s first 
noted bv lCen fc Ostrikerl <| 19991 ) and lDave etlll (I1999TI. and 
our cu rrent results are similar to those from Cen fc Ostrikerl 
(2006) using nearly the same cosmology. The addition of 
AGB feedback does not appreciably change any of these val- 



cr-derived winds have 5% fewer baryons in the WHIM 
than ^-derived winds at z = due to weaker winds at late 
times. ICen fc Ostrikerl l|2006l ) find that their baryon fraction 
in the WHIM increases by about 10% owing to galactic su- 
perwinds, which matches our results (not shown). The metal 
fraction in the WHIM is basically constant in the a-derived 
winds, in stark contrast to the "I>-derived wind model (green 
lines; also see Figure 1 of DO07) , which shows that fraction 
of metals in WHIM grow steadily from z = 4 — > 0, reaching 
25% today. The much slower cr-derived wind velocities at late 



epochs are mostly unable to shock-heat metals above tem- 
peratures where metal-line cooling dominates, and instead 
the gas cools efficiently. Although the addition of AGN feed- 
back is unlikely to change the metal content of baryons by 
more than a few percent (B07), the effect of > 1000 kms~ 
winds from QSO's such as those observed bv lTremonti et all 
(2007) could be appreciable for metals in the WHIM and hot 
phases. 

The total metal budget by baryonic phase in the second 
panel shows a minor change owing to AGB feedback, namely 
that 5% less metals are found in the galaxies, with those 
metals instead being located in the diffuse IGM. This is 
because increased star formation from gas made available 
via delayed feedback results in more winds that expel metals. 



4.2.1 Oxygen in the WHIM 

The third panel of Figure [7] shows the oxygen metallic- 
ity in various baryonic phases. Overall, oxygen metallicities 
([O/H]) remain nearly identical (within 0.1 dex) with the 
addition of AGB feedback. The 37% increase in star for- 
mation and therefore oxygen production from Type II SNe 
is counterbalanced by a 20% decrease due to the AGB pro- 
cessing of oxygen. Galactic baryons show slightly super-solar 
oxygen abundance, while the abundances in diffuse and hot 
phases are nearly one-tenth solar. 

The WHIM oxygen abundance is relatively constant 
with redshift, and shows [0/H]=— 1.76 at z = 0. Our simu- 
lations produce two distinct types of WHIM: (1) the unen- 
riched majority formed via the shock heating resulting from 
structural growth, and (2) the WHIM formed by feedback, 
which is significantly enriched. The weaker cr-derived winds 
form very little of the latter. While simulations suggest, by 
comparison with observed Ovi absorbers, that t he WHIM 
meta l licity should be around one-tenth solar (|Cen et al.l 
l200ll ; IChen et al.ll2003T ). it remains to be seen whether the 
cr-derived winds are in conflict with O VI observations. Since 
Ovi arises in both photoionized and collisionally- ionized 
gas, it could be that enough O VI is present in photoion- 
ized gas to explain the observed number density of such sys- 
tems. More over, non-equilibriu m ionization effects could be 
important (|Cen fc Fanall2006I) . A careful comparison with 
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Figure 7. The evolution of of mass and metals by baryonic 
phase in the 64/i _1 Mpc, 2 X 256 3 simulations with and with- 
out AGB feedback, and the 32/i _1 Mpc, 2 X 128 3 simulation with 
the old potential-derived winds (shown only in top two panels). 
The slower cr-derived winds at low redshift do not inject nearly as 
many metals into the WHIM. The addition of AGB feedback docs 
not change the baryonic mass fractions (Panel a) and increases 
slightly the amount of metals in galaxies (Panel b). Oxygen metal- 
licity (Panel c), which traces Type II SNe, remains nearly un- 
changed, but delayed feedback boosts carbon relative to oxygen 
(Panel d), and iron is increases at late times relative to oxygen 
(Panel e) in hot gas mostly due to Type la SNe. The global [Fe/O] 
in a test simulation without any delayed feedback is shown as the 
long dashed orange line in Panel (e). 



O VI observations is planned, but is beyond the scope of this 
paper. 

For now, we note that the oxygen abundance in the 
WHIM may be an intere s ting p robe of feedback strength. 
ISommer-Larson fc Fvnbol (|2008r ) come to the same conclu- 
sion at z = 3 when tracing oxygen content by temperature, 
noting that the stronger feedback by a top-heavy IMF pro- 
duces significantly greater amounts of oxygen in the WHIM. 
Extre mely fast winds (f wind > 1000 kms -1 ) emanating from 
AGN jTremonti et al] 120071) will also create more enriched 
WHIM. 



4.2.2 Carbon in the IGM 

In the bottom two panels of FigureQwe show the carbon and 
iron abundance relative to oxygen, a tracer of Type II SNe 
enrichment, to emphasize the differences in these two species 
influenced by delayed feedback. For carbon especially, AGB 
feedback has a significant impact. [C/O] evolution shows an 
obvious increase even at high-2, because the timescale for 
mass loss from carbon stars is 0.2 — 1 Gyr. Every phase ap- 
pears to evolve similarly with their lines sometimes blending 
in Figure [7] except the hot phase, which has at least 50% 
more carbon at z > 2 . Carbon stars lose their mass near 
sites of star formation, and this carbon is then blown out and 
shock-heated by a second generation of supernovae. The net 
effect of AGB feedback on the z = carbon metallicity is 
+0.22, +0.25, and +0.33 dex for diffuse, WHIM, and hot 
IGM respectively. This results in abundance ratios close to 
solar in all phases. 

A basic observational test of IGM enrichment models is 
the evolution of f2(Clv), i.e. the mass density in Civ sys- 
tems seen in quasar absorption lines. In Figure [5] we plot 
f2(Civ) from z — 6 — > (see OD06 for exact method of 
computing fi(Civ)) to see the effect of cr-derived winds and 
AGB feedback. In OD06 we reproduced the relative invari- 
ance in the observed trend of 0(C iv) between z = 5 — > 2 by 
counterbalancing the increasing IGM carbon content by a 
similarly lowering C IV ionization factor; the new cr-derived 
models also match the observed trend quite well, for a simi- 
lar reason. The addition of AGB feedback increases S7(Civ) 
by 70% (+0.23 dex) at z < 0.5, leading to a value consistent 
within the error bars of z s=s measurement by iFrve et ail 
(|2003l ). 

The main reason is that AGB feedback adds new fuel 
for star formation, resulting in more C IV expelled into the 
IGM at late times. Compared to the 4>-derived winds, the 
cr-derived winds push out more metals early better match - 
ing the high-z C IV observations of Ryan- We ber et al.l (|2006l ) 
an dS imcoel l|2006l ). The faster ^-derived winds at low-z raise 
the temperature of the metal-enriched IGM while pushing 
the metals to lower overdensities and lowering the C IV ion- 
ization fraction. 

The cr-derived wind model with AGB feedback is the 
first model we have explored that is able to fit the entire 
range of fl(C iv) observations from z ~ 6 — > 0. 

The cr-derived wind model with AGB feedback achieves 
higher fl(C iv) at z > 5 and z < 1 that at face value im- 
proves agreement with observations. While these data are 
uncertain and hence one should not over-interpret this im- 
proved agreement, the main point of this exercise is to show 
how our newly incorporated physical processes could have 
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Figure 8. Evolution of n(Clv) from z = 6 — » for the <3>- 
derived wind model with AGB feedback (green) and the c-derived 
wind model with and without AGB feedback (mag enta and 
blue) . Our models are compared to o bserv ations from [Songailal 
feOOll ) ( black circles'), iPettini et alj ll2003h f small filled black 
squares), iBokensberg, SargeTrt"*^Tl^au"chl ||2003|) (op en triangles), 
iFrve et alj J2003I') (open square), ISoneailal d2005l) (large open 
dots), Simco el J2006l> (large black squares), and lRvan- Weber et alj 
(2006) (black triangle is a lower limit). While all models appear 
to fit the majority of the data, there are subtle differences that 
we highlight between our models. The new tr-derived winds dis- 
tribute metals more broadly boosting O(Civ) at z > 5, while 
more carbon resulting from AGB feedback enriches the IGM in 
the local Universe. The three data points below z = 1.5 are calcu- 
lated from 8 /i -1 Mpc 2 X 128 3 simulations and from 16 /i -1 Mpc 
2 X 256 3 elsewhere, except in the case of the cr-derived winds 
without AGB feedback in which case only the small box is used. 



observational consequences. Furthermore, the results are to 
be taken with caution, owing to the small box size of our 
simulations (8 /i _1 Mpc 2 x 128 3 ) below z = 1.5; this volume 
is not nearly large enough to form the large-scale structures 
in the local Universe. It is encouraging that the values de- 
rived from this small box agree within the error bars with 
the larger U6n256vzw-(7 simulations above z = 1.5. Future 
observations at low-z by the Cosmic Origins Spectrograph, 
and at high-z with advances in near-IR spectroscopy will 
allow more relevant and detailed comparisons. 



4.2.3 Iron in the ICM 

The [Fe/O] evolution (bottom panel of Figure [7]) is domi- 
nated by Type II SNe until z ~ 2, at which point delayed 
feedback processes of Type la SNe and AGB stars become 
important. The instantaneous component of the Type la 
SNe adds 19% of the Type II SNe iron yield, which raises 
the [Fe/O] everywhere from -0.30 to -0.23 (compare to global 
[Fe/O] in a simulation without any Type la's or AGB feed- 



back . The delayed component of Type la SNe adds only 
9% more of the Type II SNe iron content generated over 
the lifetime of the Universe. However, the combination of 
the high iron yield (0.43) and the location of enrichment of- 
ten being low-metallicity regions away from the sites of star 
formation meas delayed Type la's can have a significant ob- 
servational signature in a low-density medium such as the 
ICM. The net increase of all Type la SNe (delayed and in- 
stantaneous) on the z — iron content is +0.32 dex in the 
hot component and +0.21 dex in the WHIM. AGB feedback 
increases iron content by allowing 37% more stars to form; 
this extra iron primarily remains in galaxies, but increases 
iron in the hot component by +0.15 dex and the WHIM 
by +0.07 dex. Overall, the hot iron metallicity increases by 
nearly 3x with the addition of Type la and AGB feedback. 

To demonstrate the effect on observables, we calcu- 
late the free- free emission- weighted [Fe/H] of the ICM in 
clusters/groups with temperatures in excess of 0.316 keV 
at z — in the 164n256vzw models to simulate the X- 
ray observations. We identify large bound systems in the 
simulations using a spherical overdensity algorithm (see 
iFinlator et al.ll2006l , for description). The average of over 130 
clusters/groups is -1.11 without any delayed feedback, -0.78 
with Type la SNe included, and -0.57 with AGB feedback 
also included. Hence the addition of delayed forms of feed- 
back increases the ICM [Fe/H] by 3.4x. This is now in the 
range for the canonical ICM metallicities of around 0.3 solar, 
as well as for groups between 0.316 and 3.16 keV as observed 
bv lHelsdon fc Ponmanl l|2000l ) (they found -0.60 <[Fe/H]< 
—0.36). Of course, X-ray emission at ~ 1 keV has a signifi- 
cant contribution from metal lines, and ob servations can b e 
subject to surface brightness effects (e.g. iMulchaevI .2000). 
so this comparison is only preliminary. A more thorough 
comparison of simulations to ICM X-ray observations is in 
preparation |Dave. Sivanandam fc Oppenheimerll2008l ). 

To summarize this section, we showed that Type II SNe 
remain the dominant mode of global production for each 
species we track, usually by a large margin. When consid- 
ering metals not locked up in stellar remnants (i.e. observ- 
able metals), Type la SNe only produce 22% of the cosmic 
iron and AGB stars only contribute 5% to the cosmic car- 
bon abundance. Mass feedback from long-lived stars allows 
metals to be recycled and form new generations of stars, in- 
creasing late-time star formation by ~ 30 — 40%. More im- 
portantly, the location of metals injected by delayed modes 
of feedback can significantly impact metallicities in specific 
environments. The IGM carbon content, probably the best 
current tracer of IGM metallicity, increases by 70% by z = 
when AGB feedback is included. The iron content observed 
in the ICM at least triples, primarily due to Type la SNe. 
Delayed metallicity enrichment appears to heavily affect the 
enrichment patterns of the low density gas of the IGM and 
ICM where there are relative few metals, compared to galax- 
ies where we find the metallicity signatures of Type II SNe 
dominate. Although our simulations do not produce large 
passive systems at the present epoch, it is likely that delayed 
modes of feedback will be important for setting the metal- 



2 We ran a 132nl28vzw-cr test simulation with only Type II feed- 
back, and the only major difference is the lack of iron produced 



via Type la's relative to 132nl28vzw-<r-noagb. 
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licity in and around such systems as well. Hence incorporat- 
ing delayed feedback is necessary for properly understanding 
how metals trace star formation in many well-studied envi- 
ronments. 



5 GALAXIES AND FEEDBACK 

Thus far we have examined energy balance and metallicity 
budget from a global perspective. In this section we investi- 
gate such issues from the perspective of individual galaxies. 
We will answer such questions as: What galaxies are dom- 
inating each type of feedback (mass, energy, and metallic- 
ity)? How does this evolves with redshift? Do winds actu- 
ally leave galaxy haloes and reach the IGM? What types of 
galaxies are enriching the IGM at various epochs? The key 
concept from this section is wind recycling; i.e. the prod- 
ucts of feedback do not remain in the IGM, but instead are 
either constantly cycled between the IGM and galaxies or 
never escape their parent haloes in the first place, and are 
better described as halo fountains. 

During each simulation, all particles entering a wind are 
output to a file. The originating galaxy is identified, and the 
eventual reaccretion into star-forming gas is tracked. In this 
way wind recycling can be quantified in galaxy mass and 
environment. Throughout this section we will use SKID- 
derived galaxy masses, which match our on-the-fly FOF 
galaxy finder for the vast majority of cases (cf. H2.2p . We 
use our favored cr-derived wind simulations in our following 
analysis, unless otherwise mentioned. 



5.1 Feedback as a Function of Galaxy Mass 

Figure [9] quantifies mass (upper left) , metal (upper right) 
and energy (lower panels) feedback, as a function of galaxy 
mass in our cr-derived wind simulations at four chosen red- 
shifts (2 = 6,4,2,0.5). We choose z = 0.5 to represent the 
local Universe rather than 2 — 0, because we want to fol- 
low the evolution of wind materials after they are launched 
and consider 5 Gyr a compromise as enough time for the 
winds cycle to play out, but not too much such that the 
cosmological evolution is overly significant. The upper left 
panel shows the mass loss rate in outflows as a function of 
galaxy baryonic mass. At a given galaxy mass, the outflow 
rate goes down with time, by roughly a factor of 10 from 
2 = 6^ 0.5. Remember that this is the rate of mass being 
driven from the galaxy's star-forming region; whether the 
material makes it to the IGM or remains trapped within the 
galactic halo will be examined later. 

Along with the results of our simulations, we plot two 
simple "toy models" of feedback behavior corresponding to 
galaxies forming stars at constant specific star formation 
rates (i.e. star formation rate per unit stellar mass) of 1.0 
and 0.1 Gyr -1 ; these roughly correspond to typical star 
forming galaxies at z — 6 and z — 0.5 respectively. The 
momentum-driven wind model predicts that mass feedback 



should go as M wind oc SFR/cr oc SFRxM ga | /3 if ( 
Making the reasonable assump tion that SF Roc M ga i as typi- 
cally found in simulations fe.g. |Pave]|2008l ). then this simple 
model would predict M wind oc M^H(z)~ 1/3 . The dotted 
lines in Figure [9] show these relations for our toy models. 
The red dotted line fits well to z = 6 at Af ga i = 10 9 ' 5_10 M Q 



-1/3 



showing these galaxies efficiently doubling their mass every 
1.0 Gyr, while the doubling time is around 10 Gyr for galax- 
ies at z = 0.5. 

The typical mass outflow rate reduces with time for two 
reasons: First, the star formation rates are lower owing to 
lower accretion rates from the IGM, as discussed in SH03 b 
and OD06 and as o b served by Perez-Gonzalez et al.l (|2005l ) ; 
Caput i et all (|2006l ) ; iPapovich et al.l (|2006l ). Second, galax- 
ies grow larger with time and the mass loading factors drop; 
this even despite the H{z)~ 1 ^ 3 factor that actually increases 
77 for a galaxy of the same mass at lower redshift. Hence the 
outflow rates qualitatively follow the trend seen in obser- 
vations that at high redshifts, outflows are ubiquitous and 
strong, while at the present epoch it is rare to find galaxies 
that are expelling significant amounts of mass. 

Figure [5] upper right panel, shows the mass of met- 
als launched as wind particles, which is M(Z) w i n( j oc 
M^ 3 H(z)~ 1 ^ 3 Z g3i i. For concreteness we follow the iron 
mass, although other species show similar trends; recall that 
even at z = 93% of iron is produced in Type II SNe, and 
the fraction is higher at higher redshifts. This relation can 
be thought of as the M ga i — M w i n d relation shown in the 
upper left panel convolved with the star formation- weighted 
gas mass-metallicity relations of galaxies at the chosen red- 
shifts, since it is this ga s that is being driven out in outflows. 
iFinlator fc Pavel (|2008l ) showed how our 'vzw ' model repro- 
duces the slope and scatter of t he mass-metallicit y relation- 
ship of galaxies as observed by lErb et al.l (|2006h for z = 2 
Lyman break galaxies. Again, we plot two dotted lines cor- 
responding to our toy models with aZa Ai" ga 3 dependence 
accounting for the mass-metallicity relationship normalized 
to z = 0.3 and l.0Z Q at lO n M for z = 6.0 and 0.5 re- 
spectively. The two toys models show little evolution (0.18 
dex decline from z = 6 — > 0.5, because the declining M w i n( j 
is counter-balanced by an increasing metallicity for a given 
mass galaxy toward lower redshift. In the 132n256vzw-<r sim- 
ulation, a 10 10 M Q yr _1 galaxy injects 7.7 x lO" 3 M yr~ 1 of 
iron at 2 = 6 and 2.2 x 10" 3 M o i/r~ 1 at 2 = 0.5. 

In the lower left panel we plot the total feedback en- 
ergy per time (i.e. feedback power) imparted into wind par- 
ticles as a function of galaxy mass. This power is -E w ind = 
0.5Af wind x v wind 2 oc SFRxa#(2) 1/3 oc M^H(z) 1/3 , using 
the same assumption of SFR oc M ga i- The feedback power 
is shown in units of 10 51 ergs yr -1 , which can be thought of 
as the number of SNe per year. Our simulations follow the 
trend of the toy models in terms of redshift evolution, and 
show an even tighter agreement at the low mass end versus 
Mgai than the mass feedback; the metallicity dependence 
gives greater energies to winds from lower mass, less metal- 
rich galaxies. Energy feedback is an even stronger function 
of galaxy mass than mass or metallicity feedback. 

The bottom right panel in Figure [5] shows feedback en- 
ergy relative to supernova energy (-E w i n d/-EsN). This quan- 
tity decreases with time, and increases with galaxy mass. 
The toy models shown represent Esn oc SFR oc M ga i, or 
E w i a d/Esn oc M^H(z) 1 ^ 3 . At the low-mass end, the simu- 
lations rise above the toy model and show more scatter due 
to uncertainties in the the star formation rates of the small- 
est galaxies, but the slope at most redshifts is correct. At 
the high-mass end, momentum-driven wind energy exceed 
the supernova energy, which is physically allowed since the 
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Figure 9. Feedback properties are shown as a function of galaxy mass (derived by SKID) for 4 redshifts for our cr-derived wind 
simulations at various box sizes. Red and blue dotted lines are our "toy models" calculated using momentum-driven wind relations and 
assuming specific star formation rates of 0.1 and 1.0 Gyr — , and metallicities of Z = 0.3Zq and 1.0 for z = 6 and 0.5 respectively. Mass 
and metal feedback (Panels (a) and (b)) are predicted to go as Af^ 3 , although metallicity feedback additionally depends on the gas 
mass-metallicity of galaxies, where Z oc Energy feedback, in 10 51 erg units (Panel (c)) follows the predicted Af*^ 3 relation very 



1/3 

closely, and the ratio wind energy efficiency, Ewind/^SNi (Panel (d)) rises nearly as M ^ for most redshifts. The lines on the left of 
Panel (d) show the average £? W md/^SN efficiency in the 132n256 box is declining slightly, although these values may be lower when lower 
mass galaxies not resolved in this box are included. 



UV photons produced during the entire lifetimes of massive 
stars drive winds in this scenario (see OD06, Figure 4). Sum- 
ming this ratio globally over all galaxies at each redshift, we 
obtain the values shown by the tick marks on the left side 
of the panel. Globally, the average -E W md/-EsN ratio exceeds 
unity at all redshifts (being around 1.2), and is surprisingly 
constant, declining less than 0.1 dex from z = 6 — » 0.5 in 
the new a-derived formulation of the winds. The decline of 
wind energy feedback for a given mass galaxy toward lower 
redshift is mostly counterbalanced by more massive galaxies 
driving more energetic winds at these redshifts. Less massive 
galaxies unresolved in the 132n256 box are likely to lower this 



value somewhat, so we only want to conclude that the wind 
energy is similar to the supernova energy and stays remark- 
ably unchanged with redshift in the a-derived momentum 
driven wind model. 

To summarize, the momentum-driven wind simulations 
follow trends expected from the input momentum-driven 
outflow scalings. This is of course not surprising, and at one 
level this is merely a consistency check that the new wind 
prescription and the group finder are working correctly. But 
this also gives some intuition regarding outflow properties as 
a function of galaxy mass required to achieve the successes 
enjoyed by the momentum-driven wind scenario. For exam- 
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pie, mass outflow rates should correlate with galaxy mass, 
and outflow energy in typical galaxies is comparable to, and 
perhaps exceeds, the total available supernova energy. These 
trends provide constraints on wind driving mechanisms and 
inputs to heuristic galaxy formation models such as semi- 
analytic models. 



5.2 Feedback by Volume 

We shift from examining feedback trends in individual galax- 
ies to studying feedback trends per unit volume. In order to 
facilitate observational comparisons, we use use stellar mass, 
M*, rather than baryonic mass. 

In Figure [10] we plot histograms binned in 0.1 dex in- 
tervals of the three forms of feedback at five redshift bins 
from 2 = 6^0 parameterized by the amount of feed- 
back per cubic Mpc. These histograms include all SKID- 
identified galaxies in the 8, 16, 32, and 64 h~ 1 Mpc boxes in 
order to obtain the large dynamical range covering over 5 
decades of galaxy masses. The less computationally expen- 
sive 18nl28vzw simulation was included in the histograms 
between z = — 1.5 to probe the least massive galaxies in 
this range, since the 116n256vzw run at the same resolution 
ends at z = 1.5. We also plot the median M* of a galaxy 
contributing to each type of feedback shown as vertical lines 
at the bottom of each panel. 

The mass, energy, and metal outflow rates peak at in- 
creasingly higher galaxy masses with time. The galaxy mass 
resolution limits are 1.9 x 1O 7 M at z > 4 and 1.5 x 10 s M 
below; the peaks are mostly comfortably above these resolu- 
tion limits, indicating that we have the necessary resolution 
to resolve the source of feedback across our simulation boxes 
for the three forms of feedback. The exceptions are the mass 
outflow rates at z > 2, which peak less than 1 dex from these 
limits. The changing resolution limit at z = 4 raises concern 
whether the evolution above and below this limit is real; 
however, the larger amount of evolution between z — 4 — > 1, 
despite an unchanging resolution limit shows that evolution 
at z < 4 is not just a resolution effect. 

The median galaxy M» expelling gas increases by ~ 
100 x between z = 6 — > 0, for all three forms of feedback. 
The vast majority of this growth occurs between z — 6 — > 2, 
where the median stellar mass increases by ~ 30 x in just 

2.3 Gyr. This is the epoch of peak star formation in the 
Universe, so it is not surprising that galaxies show the most 
growth in their stellar masses then. The baryonic mass (stars 
plus gas, not shown) also jumps significantly, ~ 10 — 15 x, 
but early small galaxies are more gas rich making the jump 
less extreme. Nevertheless, the evolution is much slower in 
the 10 Gyr between z = 2 — » as both median M* and M ga i 
at most triple and usually double in this longer timespan. 
This late growth could be an overestimate because of our 
overestimated star formation rates at late times in the most 
massive galaxies. It is quite possible that if our simulations 
could properly curtail massive galaxies from forming stars, 
median M* would fall toward z — 0. 

The same growth patterns for star formation-driven 
feedback are also seen in the median galaxy weight ed by 
SFR, which itself is highly correlated with M ga i as iDavel 
(2008) shows using these same simulations. Therefore the 
changing mass scales of star formation-driven feedback re- 



flect the hierarchical growth of galaxies between z = 6 ^ 0, 
with most growth occurring before z — 2. 

While it is no coincidence that all three forms of feed- 
back show similar growth rates in median M* (and M ga i), 
each feedback form has a different mass preference and peaks 
at a different epoch. Mass feedback preferentially traces 
smaller galaxies due to the inverse relationship of rj with 
Mg a i. High mass loading factors from small galaxies, which 
dominate star formation at high redshift are necessary to 
curtail the star formation density at ear ly times (OD06) and 
fit galaxy luminosity functions at z > 3 (|Finlator et al.ll2006l ; 
iBouwens etafl 120071 ) . The mass feedback density peaks at 
z ~ 3.5 at O.28M t/r _1 Mpc" 3 . The median galaxy in 
the relatively nearby Univers e bin (z = 0.5 — * 0) has 
M» = 10 10 2 M Q , or l/6th the iBell et all (|2003h M* value 
of 1O 11O M in stellar mass. 

The median M* for energy feedback is biased toward 
larger galaxies for momentum-driven winds. Our simulations 
suggest that the median galaxy adding energy to the IGM 
is an M* galaxy (M* = 1O 11 M ) in the local Universe. 
Observations disagree with this prediction, since the proto- 
typical local galaxy exhibiting feedback is more like M82, an 
~ 0.1M* galaxy. Moreover most super-M* galaxies are red 
and dead, unable to generate star formation-driven winds. 
This again suggests that a more realistic truncation of star 
formation in low-z massive galaxies could alter the exact 
values quoted here. However, at high-z where the typica l 
galaxy appears to be driving an outflow (|Erb et al.| [2006). 
these results should be more robust. The energy feedback 
density peaks at z ~ 2.5 at an equivalent energy of 5 x 10 ~ 4 
SN yr _1 Mpc -3 (10 51 ergs each). 

The metallicity feedback grows the fastest of all forms of 
feedback during the epoch of peak star formation, and does 
not obtain its maximum density until slightly past z = 1 
(M(Fe)wmd = 1-7 x lO" 4 A/ yr" 1 Mpc" 3 ). The metallicity 
feedback is simply the mass feedback modified by the mass- 
metallicity relationship of galaxies, which favors higher M t 
and lower redshift. The median M» at z = is 0.5M*, or 
~ 5x higher than the 0.1L* median galaxy as determined 
by B07. However, we do not think that this is an inconsis- 
tency due to an overestimate in massive galaxy star forma- 
tion at low redshift, but instead a result of measuring dif- 
ferent quantities. We are measuring the amount of feedback 
leaving a galaxy's star forming region whereas B07 measures 
the feedback expelled from the galaxy reaching the IGM by 
determining the effective yield as a function of v Iot . As we 
will show next, much of the mass and metals that are ex- 
pelled never leave their parent haloes and are reaccreted in 
a timescale often much less than the Hubble time. The leads 
us to introduce the important concept of wind recycling. 

5.3 Wind Recycling 

What happens to the mass and metals once they are ex- 
pelled from the galaxies' star forming regions? Is all feedback 
best described as galactic superwinds or does some feedback 
never really escape from its parent halo and should more 
accurately be considered a halo fountain? To answer such 
questions we introduce the concept of wind recycling, which 
plays an important role in feedback over cosmic time. 

By following wind particles by their particle IDs during 
the simulation run, we can track how many times the same 
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Figure 10. Mass, metallicity, and energy feedback per cubic Mpc binned by stellar mass (0.1 dex bins, derived by SKID) in 5 rcdshift 
bins covering z = 6 — > 0. Histograms include data from 8, 16, 32, and 64 /i — x Mpc in order to resolve a large range (5 decades) of galaxy 
masses. Vertical lines with colors corresponding to their redshift range show the median galaxy producing each type of feedback. For all 
forms of feedback, the median M„ increases by a factor of 100 between z = 6 — > 0; most of this increase (30 x) is between 2 = 6^2. 
However each form of feedback favors a different mass- scale: lO la2 M for mass, 10 110 Mq for energy, and 1O 1O ' 7 M for metallicity. The 
mass resolution limits are 1.9 X 1O 7 M above 2 = 4 and 1.5 X 1O 8 M below. 



SPH particle is recycled in a wind. For the 132n256vzw-cr 
simulation evolved all the way to z — we find that a wind 
particle is launched an average of 2.5 times; while 18.3% 
of SPH particles are ever launched in a wind, the summed 
number of wind launches equals 45.7% of the total num- 
ber of SPH particles. Wind recycling dominates over winds 
being launched from galaxies for the first time- the aver- 
age wind particle across time is more likely to have already 
been launched in a wind! The most important aspect is that 
wind material, once launched, cannot be assumed to be lost 
from the galaxy forever and remaining in the IGM. This is 
true despite the fact that, in our momentum-driven wind 
prescription, outflows are always ejected at speeds exceed- 
ing the escape velocity of its parent galaxy. Gravitational 
infall and hydrodynamic effects both conspire to slow down 
outflows and facilitate wind recycling. 

Figure [11] displays histograms of the number of times 
the same wind particle is recycled, iV r ec, in the 132n256 sim- 
ulation. Only 17% of all winds are particles ejected one time 
and therefore never recycled (i.e. iV roc = 0); this corresponds 
to 7.9% of all SPH particles. The record-holder is a particle 
recycled an astonishing 30 times indicating that probably in 
this case the term halo fountain may be more appropriate 
than galactic superwind. Perhaps the most telling statistic 
is half of the wind particles have been recycled 3 or more 
times. The continuous range in the number of recycling times 
blurs the distinction between a galactic superwind and a 
halo fountain, and suggests instead there is a continuum. 

The concept of recycling is not unexpected, and has 
been predicted by (|Bertone et al.ll2007h from semi-analytical 
models. DO07 showed that metals move from mean cosmic 
density at z — 2 to an overdensity of 100 by z = as baryons 
migrate into larger structures as part of cosmic structural 
growth. This means that most of the metals in the diffuse 
IGM at z ~ 2 are in galactic haloes at z ~ 0. Metals blown 
out from early galaxies to low overdensities are later reac- 
creted in the formation of larger structures, and blown out 
again. Still, the commonality of wind particles being recy- 
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Figure 11. The percentage of all SPH particles recycled iV re c 
times by 2 = in the 32h — *Mpc 2 X 256 3 simulation is shown by 
the black histogram (left scale). iV r ec = means a wind particle 
is launched only once. The red line shows the cumulative sum of 
wind particles launched per bin (right scale, normalized to 100%) 
and is the sum of (iVrec + 1) multiplied by the black histogram 
added left to right. 50% of winds are particles that have been or 
will be recycled 3 or more times indicating that wind recycling 
plays a significant if not dominant role in feedback. 



cled is surprising, especially since momentum-driven winds 
are almost always ejected at velocities well in excess of the 
escape velocity of the galaxy. It is in fact more appropri- 
ate to talk about how long it takes a wind particle to be 
recycled instead of whether it will be recycled. Metals in- 
jected by galactic superwinds cannot be assumed to remain 
permanently in the IGM; metals continuously cycle between 
galaxies and the IGM. In an upcoming paper we will quan- 



Mass, Metal, and Energy Feedback in Cosmological Simulations 21 



tify the ages of metals observed at different absorption lines 
tracing different regions of the IGM. For now we note that 
the average ages of the metals in the IGM are typically much 
shorter than the age of the Universe. 

Do wind particles generally return to their parent 
galaxy, or do they jump from galaxy to galaxy? The answer 
is the former; in the vast majority of cases a wind particle 
returns to either its parent galaxy or the result of a merger 
involving the parent galaxy. 95% of recycled wind particles 
are re-launched from a galaxy of similar or more mass than 
the previous recycling. Of course galaxies grow anyway un- 
der the hierarchical growth scenario, so a wind particle could 
join a different galaxy that has itself grown larger than its 
parent. By considering wind particles recycled within 10% of 
the Hubble time, we can diminish the bias of galaxy growth 
and explore whether winds are traveling from massive cen- 
tral galaxies to surrounding satellite galaxies possibly affect- 
ing their dynamics and enrichment histories. With this time 
limit, 97% of winds return to more massive galaxies indi- 
cating this is a less likely trend. The number jumps to 99% 
when considering winds launched from M ga i > 1O 11 M0, in- 
dicating that it is harder to escape from a more massive 
galaxy. Furthermore, a wind particle is rarely relaunched 
more than 1 comoving /i -1 Mpc apart in a 16 /i -1 Mpc test 
simulation; the minority case usually involves wind particles 
following the position of a fast moving parent galaxy. Winds 
rarely escape forever their parent halo, and winds from large 
central galaxies do not appear to disrupt satellite galaxies. 

Another way to quantify wind recycling is using the 
summed amount of material injected into winds. We intro- 
duce the quantity fi w i n d) which is the mass injected in winds, 
including recycled gas. Owing to recycling, there is no limit 
to how large this can be. Using our 132n256 run, we find 
that f2 w i n d = O.39SI5- i.e. an equivalent of 39% of the bary- 
onic mass has been blown out in a galactic superwind by 
This is significantly larger than S7, = 0.097 due to 
mass loading factors greater than one as well as recycling. 

The other cr-derived wind runs to 2 = give different 
values: f2 w ind = 0.59Qb for the 8 /i _1 Mpc box with 128 3 
SPH particles, and f2 w i nt j = 0.23f2(, for the 32 ft _1 Mpc box 
also with 128 3 particles. f2 w ind increases by 160% when mass 
resolution changes by 64 x. This jump is primarily due to a 
70% increase in 57*, but still leaves another 50% likely due 
the varying treatment of recycling at different resolutions. 
This lack of resolution convergence raises the question of 
whether the amount of recycling is set by numerics, perhaps 
arising from the poor treatment of SPH trying to model 
outflow processes occurring at resolutions below the limit of 
our simulations. 

The discontinuity in recycling rates between resolutions 
stems from how long and how far a particle reaches once it 
has been launched. Wind particles in the 32 /i -1 Mpc 2 x 256 3 
box with typical masses of 3.2 x 10 7 M Q are shot at hun- 
dreds of kms - - a huge bundle that may be more accu- 
rately thought of as a flying wind "bullet" . An actual wind 



3 Although the total number of wind launches equals 45.7% of 
all SPH particles, the typical wind particle is less massive than 
the average SPH particle, because one half of a full SPH particle 
remains when a star particle is spawned, and wind particles are 
more likely to arise from these remaining SPH particles. 



should form a bow shock resulting in a plume of material 
spreading laterally, but instead in SPH this speeding bul- 
let is slowed down viscously through the interactions with 
neighboring particles. The number of interacting neighbors 
over the same physical distance is lower at less resolution, 
and a wind particle will take longer to slow down; this is 
an unfortunate but unavoidable consequence of SPH. Hence 
lower resolution simulations inhibit recycling as evidenced 
by the average number of times an individual wind par- 
ticle is launched: 3.0, 2.5, and 2.1 times for 4.7, 38, and 
300xl0 6 M Q SPH particle resolution (to 2 = 0). This lack of 
resolution convergence is quantifiable by the median recy- 
cling timescale of wind particles in different resolution sim- 
ulations as shown in the upper right panel of Figure [121 and 
is discussed in the following subsection. Increased resolution 
suggests that wind recycling should move f2 w i n d to an even 
higher value, probably exceeding 1/2$1{,. Given that galaxies 
below a certain mass cannot form owing to the presence of 
an ionizing background, there will be a limit to how high 
f^wind can be; in the future we hope to run simulations that 
can achieve such resolutions in a cosmologically representa- 
tive volume. For now, the wind recycling predictions should 
be considered as illustrative rather than quantitative. 

5.4 Wind Recycling Timescales 

We now examine wind recycling as a function of galaxy 
mass. Specifically, we want to know what galaxies are able 
to inject their metals into the IGM. The common thought 
is that low-mass galaxies lose their metals more easily be- 
cause winds ca n escape from these gala xies' shallower poten- 
tial wells (e.g. iDekel fc SilkHl986r i. As lDekel fe Wool (|2003h 
showed, if the outflow energy couples efficiently to ambi- 
ent halo gas, this can yield the mass-metallicity relation 
as observ ed, as well as other pro perties of dwarf galaxies. 
However, iFinlator fc Pavel ([2008 ) showed that a model in 
which galaxies expel material at constant velocity (i.e. the 
"constant wind" model of OD06) does not reproduce the 
observed mass-metallicity relation, primarily because out- 
flows do not in fact couple their energy efficiently to ambi- 
ent gas. Instead, the observed mass-m etallicity relation (e.g. 
lErb et alJlbood : iTremonti et al.ll2004r i are better reproduced 
in the momentum-driven wind scenario. In this case, since 
"wind oc y/— Sgai (roughly), all outflows have an approx- 
imately equal probability of escaping their parent haloes 
independent of mass, when considering only gravitational 
interactions with an isolated halo. 

The left panel of Figure [12] shows the virial ratio, 
E T!i t = 0.5u w ind 2 /$, as a function of galaxy baryonic mass, 
for winds launched at four redshifts. At any given redshift, 
the ratio is roughly constant, but rises slightly at low masses 
and then declines to higher masses. At low masses, the trend 
arises from the mild metallicity dependence of v w ind together 
with the mass-metallicity relation. Moving to higher masses, 
since the potential $ includes the galaxy potential and the 
environmental potential of the group/cluster in which it lies, 
and since massive galaxies live in denser environments, $ 
is greater. At the highest masses, our imposed speed limit 
reflecting the assumption that wind energy cannot exceed 
twice the supernova energy reduces wind velocities. Note 
that by definition, <E>-derived winds will have a constant _E ra t 
for a given launch redshift (modulo wind speed limits); hence 



22 B. D. Oppenheimer & R. Dave 
1.5 



x 
in 
d 

I 

O 



0.5 



-0.5 




10 

log(M gal ) 



10 



9.5 



CO 



o 



8.5 



~i 1 r 



n 1 r 



(b) 




l I l l l 



J_ 



8 10 12 

log(M gal ) 



Figure 12. The virial ratio (-E ra t, the ratio of wind kinetic energy to potential well energy) are displayed as a function of galaxy mass 
(derived by SKID) at 4 launch redshifts in Panel (a) for our cr-derived wind simulations at various box sizes. The shallow turnover at 
high-mass is mostly due to imposed wind speed limits. Panel (b) shows that the recycling timescale (the median time for a wind particle to 
be relaunched in a wind) is a much stronger function of M ga i. After exploring alternatives for the M ga i — t rGC anti-correlation, we conclude 
winds from more massive galaxies encounter a denser environment slowing the wind particle faster and allowing it to be re-accreted faster 
The two dotted lines show the galaxy mass dependence of t 
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is dominated by slowing due to environment, M ^ ' , or the gravitational 



—1/2 

free-fall time, M ^ ; wind particles spend most of their time falling back into a galaxy. The poor resolution convergence in t 



rec indicates 

lower resolution simulations overpredict the t Iec and underpredict the total amount of wind recycling. To see wind recycling in action in 
our simulations, see the movies and an explanation at http://luca.as.arizona.edu/~oppen/IGM/rccycling.html 



the curvature towards high masses is predominantly a reflec- 
tion of cr-derived winds. 

In practice, outflows are not only confined by the grav- 
itational potential of the host galaxy; cosmic infall and hy- 
drodyn amic effects can als o be quite important. For in- 
stance, iFerrara et all (|2005t ) used a simulation of a typical 
LBG at z = 3 (M dyn = 2 x lO n M ) to show that metals 
never reach the IGM and instead are confined to the sur- 
rounding hot halo gas by infalling gas that creates a shock 
interface with the outflow. 

One way to quantify such additional effects is to con- 
sider a new timescale called the recycling timescale, t rcc , 
which is defined as the median time for a wind particle 
to be re-ejected again as a wind particle, or else fully con- 
verted into a star. At all launch redshifts explored, for all 
Mg a i , more than half of wind particles recycle meaning that 
t ICC can be measured. The right panel of Figure [12] shows 
tree — Mg a i relation, demonstrating a strong anti-correlation 
of ircc with galaxy mass. The simple fact that winds recycle 
on short timescales even though they are launched at ve- 
locities well above the escape velocity is an indication that 
slowing by ambient gas and surrounding potentials is signif- 
icant- 

If the galaxy's gravity is dominant in confining its out- 
flow, one would expect that the time spent away from 
the galaxy would be approximated by twice the free-fall 
timescale, iff (one ia outward and another one back in, or 



For new movies showing wind recycling 
in action in our simulations, please visit 
http://luca.as.arizona.edu/~oppen/IGM/recycling.html 



simply the orbital timescale for an orbit of eccentricity 1). 
ts scales as Rtiru^dyl 2 ' where i?t U m is the turn-around dis- 
tance of the wind particle, and Md yn is the dynamical mass 
of the galaxy treated as a point source for the sake of simplic- 
ity. Following Newtonian dynamics, -Rtum oc A'/dyn/">wind~, 

1/3 

and the i> W md oc M ^ relation for momentum-driven winds, 
remains invariant as a function of galaxy mass while 
Rtum oc M 1//3 . Momentum-driven winds should reach a 
maximum distance from their parent galaxy that is a con- 
stant multiple of the virial radius since u w ind scales with the 
virial energy of an isolated halo. The strong trend of tree 
with Mgai indicates that larger-scale potentials and hydro- 
dynamic effects are dominant. 

Is the particle spending more time in the IGM/galactic 
halo or in the galactic ISM before being relaunched? The an- 
swer is that the vast majority of the time between recycling 
is spent outside the star forming regions of a galaxy. This 
is not surprising, because in our simulations the star forma- 
tion timescale is tied completely to the accretion timescale; 
once a particle is in a galaxy it gets converted to a star 
or blown out in a wind much more quickly than it was 
accreted |Finlator fc Davil200ct ). Therefore, we subdivide 
t roc into two timescales, one leaving the galaxy to reach 
Rtum while being slowed by hydrodynamical forces, tout 
and the other returning to the galaxy in approximately one 
iff. For simplicity, let's say the hydrodynamic forces slow 
the wind particle down at a constant rate, in which case 
tout = .Rturn/(«wmd/2). The total recycling timescale can be 
approximated as 

tree f Rtum \ f 100 kmS' 1 ^ 

lGyr ■ ^ lOOkpc J \ w wind J 
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if we assume that Md yn = !2m /Qb X M ga i and J?tm-n is in 
physical kpc. We can now solve for 7?tum by taking the tree 
from Figure [12] and the average v W md from Figure [3] i?turn is 
just an upper limit for how far wind particles with median 
t rcc can extend from a parent galaxy, because in reality a 
particle does not likely go out and then immediately fall 
back into a galaxy in two steps; wind particles often spend 
time orbiting around their parent galaxies as well. 

In principle, 7?tum can be tracked directly in simula- 
tions, however our current suite of runs did not output posi- 
tion information of wind particles owing to the large storage 
requirements. We have checked for isolated cases that our 
simple formula (eq. [9} yields roughly correct results. 

We plot i?turn in the left of Figure 1131 focusing first 
on the trend over 2.5 decades of galaxy masses at z — 0.5, 
which we consider here the local Universe since trends do not 
evolve much until z = 0. Rtum stays nearly constant from 
the size scale of dwarf galaxies to galaxies above L* while 
t re c declines. Dotted lines correspond to the radius, r2oo, of 
a NFW halo (|Navarro. Frenk. fe Whltl 19971 ). The point at 
which J?tum and V2oo intersect is the approximate transition 
mass below which winds reach the IGM and above which 
winds are confined to their haloes. Metals rarely enter the 
IGM from galaxies above 10 n MQ in the local Universe. Sev- 
eral significant conclusions can be drawn from this behavior 
at low- z. 

First let us consider our wind model and the Milky Way, 
which has M ga i ~ 10 n A/Q. Our wind model produces an 
outflow from a Milky Way-type galaxy, with initial speeds 
of a few hundred kms -1 . One might claim that this imme- 
diately invalidates our model, since the Milky Way is not 
observed to have an outflow. Yet quantitatively, our model 
predicts that a Milky Way-sized galaxy should have a recy- 
cling time of less than 1 Gyr, and the typical wind particle 
will not venture beyond 85 kpc (if we assume not much 
has changed between z = 0.5 — > 0). In other words, in our 
model, the Milky Way is not driving a classical outflow as 
seen from local starbursts, but rather is continually sending 
material up into the halo and having it rain back down on 
a timescale of 1 Gyr. We call this a halo fountain, in anal- 
ogy with a galactic fountain that operates on smaller scales. 
Indeed, we speculate that there may not be a fundamen- 
tal difference between galactic fountains and halo fountains, 
but rather gas is being thrown out of the disk at a range 
of velocities; however, our simulations lack the resolution to 
address this issue directly. 

How might such a halo fountain be observed? One possi- 
bility is t hat it has already been se e n, as high velocity clouds 
fflVC's). lw"akker fc van Woerdenl (|l997f) predict that given 
the observed rate of material going into HVC's of about 5 
M02/r _1 then it should take ~ 1 Gyr for all the gas in the 
ISM to cycle through this halo fountain. Hence the Galaxy 
may have an active halo fountain recycling its material on 
a timescale much less than a Hubble time, despite not re- 
sembling anything like a starburst galaxy. The difficulty of 
observation of feedback within our own Galaxy may just in- 
dicate how invisible yet ubiquitous galactic-scale winds are. 

Secondly, with Rt uln nearly constant, t Tec should be pro- 
portional to M~^ 3 if dominated by t out and M~^ 2 if tg is 



larger. The black dotted line in the right of Figure [12] corre- 
sponding to the latter case (the steeper of the two) appears 
to more closely match the general trends indicating that 
tff dominates; wind particles spend a majority of t rcc falling 
into galaxies (Panel (b) in Figure [T3"1l . Our calculations above 
agree grows larger than t out for smaller galaxies, but the 
two timescales are similar for extremely massive galaxies, 
Mgai ~ 1O 12 M . 

Thirdly, smaller galaxies live in less dense environments 
where hydrodynamic slowing takes longer. To demonstrate 
this we look at another parameter, the minimum overden- 
sity reached by wind particles, which should approximately 
correspond to the density at -Rtum- For the subset of wind 
particles we track the minimum density achieved before re- 
cycling, p m in. Figure PHI (solid lines) shows the median p m in 
as a function of the baryonic mass of the originating galaxy, 
Mg a i, from the 132n256 run. At all redshifts, smaller galax- 
ies push their winds to lower densities, consistent with them 
having longer recycling times. Their less dense environments 
slow winds over a longer tout , an d allow them to reach a sim- 
ilar i?turn as more massive galaxies despite lower « w i n d- We 
also plot a long-dashed line that show the average density 
within 1 comoving Mpc sphere around the z — 0.5 galaxies. 
These show a similar trend, indicating that environmental 
dependence is the primary factor in how far winds reach 
into the IGM and how long they remain there. The fact 
that p m in is much higher than the density within 1 Mpc is 
an indication that these winds are traveling much less than 
that distance as our calculations above indicate (~ 80 kpc 
comoving at z = 0.5). 

Although a possible explanation for the trend in p m i n 
is larger haloes have higher densities at the same distance, 
this cannot account for the dependence. The density of NFW 
haloes decline nearly as 1/r beyond r2oo resulting in a p m i n oc 
Mgai dependence much steeper than plotted in Figure 1141 
below z — 2; for all but the most massive galaxies the density 
contribution of the parent halo is much smaller at p m i n . The 
flattening of this relation toward z — at M ga i < lO n M0 
means these galaxies more likely live in denser environments. 
For more massive galaxies, the parent halo itself is more 
responsible for the hydrodynamical slowing. 

Finally, small galaxies can enrich a similar volume as 
large galaxies leading to major implications as to which 
galaxies enrich the IGM. Metals from small galaxies have the 
advantage of staying in the IGM longer, however the disad- 
vantage of their winds being less metal enriched (e.g. Panel 
(b) of Figure [9]). We plan to look at the origin of IGM metal 
absorbers in a future paper, exploring what galaxies enrich 
the IGM, how long these metals have been in the IGM, and 
what distance metals travel from their parent galaxy. 

As a side note, we admit that the global SFRD is over- 
estimated at late times (z < 0.5); the 132n256 run show- 
ing nearly 4x as much SF at z = relative to observa- 
tions (Figure 2| . The primary reason for this discrepancy is 
the lack of a quenching mechanism in the largest galaxies; 
galaxies above 1O 11 M0 account for 54% of star formation 
at z = 0, and their continued late growth creates far too 
many galaxies above this ma ss as compared to observations 
l|Perez-Gonzalez et al.l 120081 ) . The influence of these over- 
massive galaxies on the IGM however appears to be remote 
due to their inability to inject metals beyond their haloes 
and their short recycling timescales. 
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Figure 13. Panel (a) shows Hturn, the maximum radial extent a median wind particle may extend from its parent galaxy calculated from 
equation[9]in the 132n256 simulation. -Rturn generally is larger at smaller M ga i; this trend coupled with the fact t hat t rec is longer indicates 
that s maller galaxies can more easily enrich the IGM. Dotted lines corresponding to a NFW halo radius, 7*200 iNavarro. Frenk, fc White] 
l|l997h . show that winds from L* galaxies do not escape their parent haloes at z < 1. Panel (b) shows that winds from small galaxies 
spend more time on their journey returning to a galaxy than being blown out. 
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Figure 14. The minimum density a wind particle achieves be- 
tween recyclings is shown to be a strong function of galaxy mass 
at all redshifts (colored lines). Lower mass galaxies have longer 
recycling times, because winds are launched into less dense en- 
vironments where shocks take longer to slow and turn around 
winds. The environmental density within sphere of 1 comoving 
Mpc as a function of M ga i at z = 0.5 (black dashed line from 
132n256 simulation) shows the same trend as p m i n /p indicating 
that environment is the primary factor in how far winds reach 
into the IGM. The NFW halo profile alone cannot explain this 
trend as halo profiles drop much more sharply and the environ- 
mental density dominates for all but the most massive galaxies 
at low- z. 



Turning to the evolution of recycling, t lec grows moder- 
ately longer at lower redshift for a given M g ai. For a 10 10 M Q 
galaxy, t ICC is 1.1, 1.5, 1.6, and 2.0 Gyr for launch redshifts 6, 
4, 2, and 0.5 respectively in the 32 /i _1 Mpc box. This trend is 
in place despite _E rat declining by 1 decade from z — 6 — > 0.5 
(see Figure [T27 a)); wind kinetic energy is declining relative 
to the potential at their launch location. Again, the overrid- 
ing variable is the slowing of the wind by the environment; 
the average physical density declines by a factor of 100 in 
this interval, making it easier for winds to travel further 
despite a factor of 10 x less energy input into the winds. 
Galaxies of the same mass are more likely to live at higher 
overdensities at lower redshift as structural growth makes 
more groups and clusters, but the declining physical den- 
sities and slower Hubble expansion of the Universe toward 
low redshift outweigh this. 

While the t lec grows at low redshift, the amount of the 
Universe it enriches sharply declines. The calculations of 
fltum for equation [9] shows a moderately increasing phys- 
ical distance from z = 6 — > 0.5 with values of 56, 71, 75, and 
80 kpc for the four respective launch redshifts when con- 
sidering the same 1O 1O A40 galaxy. However, the comoving 
volume such a galaxy is enriching is 35x greater at z = 6 
than at z = 0.5. Momentum-driven winds can more eas- 
ily expel metals into the IGM at high redshift, whereas the 
overdense environments combined with weaker wind veloci- 
ties and the overall larger scale of the Universe leave metals 
perpetually recycling in halo fountains in groups and clus- 
ters at low redshift. These are the primary reasons metals 
migrate from the IGM to galaxies as the Universe evolves. 

Our results compare favorably to those of lBertone et al.1 

l|2007h . who follow energy-driven winds (77 oc a 2 and 
«wind oc a) in their semi-analytical implementation using a 
two-phase process for the d ynamical wind evolution (see also 
|Pe Lucia fc BlaizorJ l|2007h ). At z = they also find that it 
is the smaller galaxies (10 10 - 5 Af Q < M dyn < 10 11,5 M Q , i.e. 
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~ 1/10M*) that most efficiently inject their metals into the 
IGM (~ 40% of the time), while metals from more massive 
haloes than 1O 12 M0 rarely escape the halo and are recycled 
(see their Figure 9). This is despite their u w i n( j-dependence 
(same as ours but without the scatter) exceeding the es- 
cape velocity. The redshift dependence is harder to compare 
to our results due to differing definitions of escaped frac- 
tions into the IGM, but they do also see smaller haloes, the 
progenitors to massive galaxies today, as the most efficient 
enrichers of the IGM up to z = 4.5. 

Lastly, t roc and p m i n show similar resolution convergence 
issues as fi W md, for similar reasons. Lines at the same red- 
shift do not overlap in either of these plots, indicating the 
t roc is longer for the same M ga i at lower resolution. As dis- 
cussed before, coarser resolution simulations slows down par- 
ticles over a longer physical distance allowing wind particles 
to travel artificially too far at low resolution; this is why 
there is less recycling at lower resolution. For well-resolved 
galaxies, the lines of different resolution appear to nearly 
converge. The hydrodynamical treatment of wind particles 
after they are launched is not nearly as resolution-converged 
as the feedback properties on per galaxy basis (see Figure 
[9|. Recycling times are overestimated when under-resolved. 
Hence recycling times may actually be shorter than we've 
shown. 

In summary, we have shown that wind recycling is an 
important phenomenon for understanding the evolution of 
galaxies and the IGM. Winds from galaxies are typically ex- 
pelled and re-accreted on timescales short compared to a 
Hubble time. This occurs despite the fact that winds are 
typically launched with plenty of kinetic energy to escape 
its halo. Such gravitational arguments are therefore not very 
relevant to understand how material cycles through galax- 
ies and the IGM; instead, infall and shocks generated by 
outflows are more important in setting the wind distribu- 
tion length- and time-scales. The recycling time is there- 
fore strongly anti-correlated with galaxy mass, owing to the 
fact that more massive galaxies live in denser environments. 
These denser environments make it harder for winds to 
travel beyond their parent halo at late times causing the mi- 
gration of metals from the IGM to galaxies from z = 2 — > 0. 
At the present epoch, small galaxies can still expel material, 
but larger galaxies are more aptly described as having halo 
fountains, in which material is constantly kicked up into the 
halo before raining down. 



6 SUMMARY 

We introduce a new version of Gadget-2, with improve- 
ments designed to explore mass, metal, and energy feed- 
back from galactic outflows across all cosmic epochs. We 
add two major modules designed to make the code better 
suited to explore the low-z Universe: (1) A sophisticated 
enrichment model tracking four elements individually from 
Type II SNe, Type la SNe, and AGB stars; and (2) an on- 
the-fly galaxy finder used to derive momentum-driven wind 
parameters based on a galaxy's mass. 

We first run test simulations to explore global energy 
and enrichment properties with and without AGB and Type 
la feedback, and with our old and new wind implemen- 
tations. Focusing on our new (galaxy mass-derived) winds 



including all sources of feedback, we then run several 34 
million-particle simulations to explore feedback over the his- 
tory of the Universe and over a large dynamic range in 
galaxy mass. We also track a representative subset of wind 
particles to study in detail how mass, metals, and energy 
are distributed by outflows. 

Our new chemical enrichment model enables us to in- 
vestigate the global production and distribution of key indi- 
vidual metal species. Globally, metal production of all four 
species tracked (C, O, Si, Fe) is dominated by Type II SNe at 
all redshifts. Type la SNe add significantly to the iron con- 
tent of hot, intracluster gas, especially at z < 1. AGB stars 
add moderately to the IGM carbon abundance by z — 0, and 
provide fresh (enriched) gas for recycling into stars which in- 
creases global star formation at later epochs. Carbon yields 
from AGB stars cannot be ignored even at high-z, because 
carbon AGB stars enrich on timescales much less than a 
local Hubble time (i.e. 200 Myr-1 Gyr). Due to the com- 
plex interplay between instantaneous and delayed recycling 
from various forms, metallicity patterns in the IGM and 
ICM cannot be straightforwardly used to infer the enrich- 
ment patterns in the host galaxies responsible for polluting 
intergalactic gas. 

We study enrichment patterns subdivided by baryonic 
phase. The total metal mass density is split roughly equally 
between galaxies and the IGM from z = 6 — > 2. Shocked 
intergalactic gas (WHIM and ICM) contains a fairly small 
portion of the global metal mass at all epochs. At z < 2, 
metals tend to migrate from the IGM into galaxies, so that 
by z — about two-thirds of the metals are in galaxies 
(i.e. stars and cold gas). The combination of increased car- 
bon and gas recycling from AGB feedback results in the 
IGM being significantly more carbon-enriched, which helps 
reproduce the relatively high observed mass density in C IV 
absorption systems at z ~ 0. 

Our new galaxy mass-based wind model im- 
plementation provides a more faithful represen- 
tation of the momentum- d riven wind model of 
iMurrav. Quatert. fe Thompson! ([2005) , and yields out- 
flows that are in better broad agreement with observed 
outflows. In particular, our implementation results in faster 
winds at high-z and slower wind speeds at low-z compared 
to our old local potential-derived winds. The fast early 
winds are able to enrich the IGM at early times as observed, 
while the slow late winds mean that most galaxies today 
are not driving material into the IGM at all. Qualitatively, 
this better agrees with observations indi cating that mos t 
galaxies at z ~ 2 — 3 drive powerful winds l|Erb et al.ll2006l ) , 
while today galaxies are rarely seen to have strong outflows. 

We examine bulk properties of outflows as a function 
of galaxy mass. We find that mass, metallicity, and en- 
ergy feedback as a function of galaxy baryonic mass roughly 
follow the trends predicted by momentum-driven winds: 
Mwmd oc M^H-^iz) and £ wind oc M^H x ' 3 {z). The 
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metal outflow rate is M w i n dZ g ai, where Z ga i oc M ^ with a 
proportionality constant that i ncreases with time as giv en by 
mass-metallicity evolution (cf. iFinlator fe D ave 2008) . The 
stellar mass of the typical (median) galaxy most responsi- 
ble for each particular form of feedback increases by a fac- 
tor of ~ 30 between z — 6 — + 2, but only ~ x2 — 3 from 
z — 2 — > 0. Each form of feedback traces a different mass 
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scale: mass feedback ~ \M* , energy feedback ~ M* , and 
metallicity feed back ~ ^M*. Mea suring these characteris- 
tic masses (e.g. iBouche et aDl2007l ) offers the possibility to 
test whether momentum-driven wind scalings are followed 
globally. 

The wind energy relative to the supernova energy scales 
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roughly as M'^ , and is within a factor of a few of unity 
at all epochs and galaxy masses. Given expected radiative 
losses from SN heat input into the ISM, it seems that SNe 
will have a difficult time providing enough energy to pol- 
lute the Universe as observed. An alternate source of energy 
would ease this tension, such as photons from young stars 
whose total energy can exceed the supernova energy by sev- 
eral orders of magnitude. This adds to the circumstantial 
evidence supporting the idea that galactic outflows may be 
driven in large part by radiation pressure, as postulated in 
the momentum-driven wind model. 

We find that wind recycling, material ejected as outflows 
and then re-accreted and ejected again, turns out to be a 
remarkably common occurrence with significant dynamical 
repercussions. By following individual wind particles in our 
simulations down to z = 0, we find that multiple recyclings 
are the norm, and that the typical wind particle has been 
ejected three to four times. In other words, outflow mate- 
rial being reaccreted onto a galaxy dominates over outflow 
material launched into the IGM forever. Since in our wind 
model all outflows are launched at speeds well above the 
escape velocity of galaxies, this indicates that outflows are 
mainly slowed through hydrodynamic interactions, allowing 
them to rejoin the hierarchical accretion flow into galaxies. 
Approximately 20% of the baryons participate in an out- 
flow, but owing to multiple launchings the aggregate mass 
of baryons ejected exceeds half of the total baryonic mass. 
The two key corollaries of wind recycling as seen in our sim- 
ulations are: (1) Material driven in an outflow cannot be 
assumed to remain in the IGM forever, and (2) Gravita- 
tional energetic considerations are generally not relevant for 
determining how far outflow propagate into the IGM. 

We examine wind recycling as a function of galaxy mass. 
The recycling time scales roughly as t rcc oc Af_j , as ex- 
pected if environmental effects are dominating the retarda- 
tion of outflows. Larger galaxies have shorter recycling times 
despite launching winds at larger speeds, because they live in 
denser environments; the minimum overdensity achieved by 
winds scales with galaxy mass. As expected from the ubiq- 
uity of multiple recyclings, the recycling time is generally 
fairly short, roughly io 9±0,5 years, increasing only mildly at 
lower redshift. Our analysis suggest that winds generally re- 
turn to the galaxy from which they were launched (or its 
descendant), though a full merger tree construction is re- 
quired to confirm this. 

It is possible to estimate how far outflows travel from 
their host galaxies (-Rturn), as a function of galaxy mass and 
redshift. Remarkably, Rtum = 80±20 physical kpc at all red- 
shifts and masses. There are weak trends for higher Rtum at 
smaller masses and lower redshifts; both are consistent with 
ambient density being a key determinant for how far winds 
travel. The constant physical distance means that outflows 
at early epochs are able to enrich a significant fraction of the 
Universe, while outflows at later epochs are more confined 
around galaxies. This is the reason metals migrate from the 
IGM to galaxies between z — 2 — > 0. 



Comparing -Rtum to halo radii, we see that at z > 2, typ- 
ical L*-sized galaxies have outflows that escape their host 
halos into the IGM, while at z < 1 outflows are generally 
confined within galactic halos. This gives rise to the concept 
of halo fountains, where low-z galaxies are constantly kick- 
ing gas out of their ISM into the halo but no further, and this 
material rains back down onto the ISM on timescales of order 
1 Gyr or less. If correct, then even galaxies not canonically 
identified as having outflows (such as the Milky Way) may 
in fact be moving a significant amount of material around its 
halo. This halo fountain gas would be quite difficult to de- 
tect, as it is likely to be tenuous and multi-phase; we broadly 
speculate it might be resp onsible for high-velocit y clouds or 
halo Mgll absorbers (e.g. iKacprzak et alj |200S| ) . We leave 
a more thorough investigation of the observational conse- 
quences of halo fountains for the future. 

As a final caveat, it should be pointed out that detailed 
properties of how outflows propagate out of galaxies are not 
as well-converged with numerical resolution as we would like. 
This may be due to our particular way of implementing out- 
flows in a Monte Carlo fashion combined with difficulties of 
SPH in handling individual outflowing particles. Here we 
have focused on qualitative trends that appear to be robust 
within our limited exploration of resolution convergence, 
without making overly detailed quantitative predictions. In 
the future we plan to investigate how to implement outflows 
in a more robust way within the framework of cosmological 
hydrodynamic simulations. 

Wind recycling and halo fountains provide two new 
twists on the idea of galactic outflows. If our models are cor- 
rect, then there is a continuum of outflow properties from 
galactic fountains that barely kick gas out of the disk, to 
halo fountains that cycle material through gaseous halos of 
galaxies, to large-scale outflows that are e.g. responsible for 
enriching the IGM. In short, outflows may be considerably 
more ubiquitous and complicated than previously thought, 
and hence understanding their effects on galaxies as a func- 
tion of mass, environment, and epoch will be even more crit- 
ical for developing a comprehensive model for how galaxies 
form and evolve. 
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